STRESS  ANALYSIS  OF  A  DOUBLE-DOUBLER  JOINT 
IN  COMPOSITE  STRUCTURE 


By 

MYUNG  KYUN  PARK 


A  DISSERTATION  PRESENTED  TO  THE  GRADUATE  SCHOOL 
OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 
DOCTOR  OF  PHILOSOPHY 

UNIVERSITY  OF  FLORIDA 


1990 


ACKNOWLEDGEMENTS 

The  author  deeply  appreciates  the  efforts  of  his 
committee  chairman,  Professor  L.  E.  Malvern.  For  his  advice, 
guidance,  and  encouragement,  the  author  is  greatly  indebted. 

The  author  is  grateful  to  Professor  B.  V.  Sankar 
especially  for  suggesting  study  of  the  stress  singularity 
treatment  and  for  helpful  suggestions  in  carrying  out  all  the 
analyses. 

Special  appreciation  is  extended  to  Professors  C.  T.  Sun, 
S.  Y.  Lu,  and  Professor  David  C.  Wilson  for  their  assistance 
and  serving  on  the  supervisory  committee. 

The  author  thanks  the  Department  of  Aerospace 
Engineering,  Mechanics  and  Engineering  Science  and  the  Center 
for  Studies  of  Advanced  Structural  Composites  of  the 
University  of  Florida  for  their  assistance  and  financial 
support  on  this  research. 

The  author  thanks  all  his  family  members  for  their 
encouragement  on  his  pursuing  an  advanced  degree  in  the  United 
States. 

Finally,  true  thanks  are  due  to  the  author's  wife,  Boeun, 
who  with  her  prayers,  love  and  forbearance  supported  the 
author  to  finish  his  study  at  the  University  of  Florida. 

ii 


TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGEMENTS    ii 

LIST  OF  TABLES    vi 

LIST  OF  FIGURES    vii 

ABSTRACT    XV 

CHAPTERS 

1  INTRODUCTION    1 

1 . 1  Introduction    1 

1.2  Objective  and  Scope    2 

2  REVIEW  OF  PREVIOUS  WORK   5 

2 . 1  Introduction    5 

2 . 2  Analytical  Study    7 

2.2.1  Volkersen ' s  Theory   "  ]  7 

2.2.2  Goland  and  Reissner  Theory    9 

2.2.3  Plantema's  Analysis   13 

2.2.4  Kelsey  and  Benson  Theory    14 

2.2.5  Pahoja's  Analysis    15 

2.2.6  Renton  and  Vinson's  Analysis    19 

2.2.7  Hart-Smith's  Analysis    21 

2.3  Finite  Element  Studies    23 

2.3.1  Introduction   \  23 

2.3.2  Ahluwalia's  Analysis  .]  24 

2.3.3  Wooley  and  Carver's  Analysis    25 

2.4  Experimental  Studies    26 

2.4.1  Mclaren  and  Maclnnes's 
Investigation    26 

2.4.2  Tuzi  and  Shimada's 

Investigation    27 

3  LINEAR  ELASTIC  STRESS  ANALYSIS    29 

3.1    Introduction    29 

3  . 2    Joint  Description    30 

3  .  3  Determining  Material  Properties  ." .' ." ." 33 
3 . 4    Finite  Element  Model  Representation ' ." .' ."  .* .' 36 

iii 


3 . 5    Finite  Element  Results 


40 


3.5.1  Case  1 .  Aluminum  Adherends  and 
Epoxy  adhesive  in 

Simple  Tension    43 

3.5.2  Case  2.  Aluminum  Central 
Adherend,  Epoxy  Adhesive,  and 
Composite  Outer  Adherend 

in  Simple  Tension    52 

3.5.3  Case  3 .  Aluminum  Adherends  and 

Epoxy  Adhesive  in  Pure  Bending   53 

3.5.4  Case  4 .  Aluminum  Central 
Adherends,  Epoxy  Adhesive,  and 
Composite  Outer  Adherend 

in  Pure  Bending   68 

4  PARAMETRIC  STUDY  OF  THE 

DOUBLE-DOUBLER  JOINT    78 

4.1  Introduction    78 

4.2  The  Material  Parameter   79 

4.2.1    Influence  of  E/Ea   80 

4.3  The  Geometric  Parameters    85 

4.3.1    Influence  of  Overlap  Length    87 

5  TREATMENT  OF  STRESS  SINGULARITY 
AT  BONDING  CORNERS  IN  THE  DOUBLE- 
DOUBLER  JOINT    89 

5.1  Introduction    89 

5.2  Stress  Singularity  at  Interface 

Corners  in  Bonded  Dissimilar  Materials  ...  90 

5.3  Calculation  of  Order  and  Intensities 
of  Stress  Singularities  at  the  Bonding 
Corners  of  a  Double-doubler  Joint    98 

5.3.1  Case  1.  Half  Plane  Bonded 

to  a  Quarter  Plane    98 

5.3.2  Case  2.  Three  Quarter  Plane 

Bonded  to  a  Quarter  Plane    103 

5.4  Specified  Boundary  Displacement  Method  ...  106 

5.4.1  Introduction    ^06 

5.4.2  Procedure   [\  107 

5.4.3  Benefits  of  Submodelling  no 

5.4.4  Coarse  Model  and  Submodel   no 

6  RESULTS  AND  DISCUSSION    117 

6.1  Introduction    117 

6 . 2  Summary  and  Conclusions   I '.  ns 

6.3  Recommendations  for  Future  Work  121 

iv 


APPENDICES 


A     CASE  1.  THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  DOUBLE-DOUBLER  REINFORCEMENT  OF  A 
CONTINOUS  MEMBER  AT  DIFFERENT 

LEVELS  OF  Y   124 

B     CASE  1.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  A  DOUBLE-DOUBLER  JOINT  AT 

DIFFERENT  LEVELS  OF  Y    129 

C     CASE  2.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  DOUBLE-DOUBLER  REINFORCEMENT  OF  A 
CONTINOUS  MEMBER  AT  DIFFERENT 

LEVELS  OF  Y    134 

D     CASE  2.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  A  DOUBLE-DOUBLER  JOINT  AT 

DIFFERENT  LEVELS  OF  Y    13 9 

E     CASE  3.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  DOUBLE-DOUBLER  REINFORCEMENT  OF  A 
CONTINOUS  MEMBER  AT  DIFFERENT 

LEVELS  OF  Y    144 

F     CASE  3.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  A  DOUBLE-DOUBLER  JOINT  AT 

DIFFERENT  LEVELS  OF  Y    15 3 

G     CASE  4.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  DOUBLE-DOUBLER  REINFORCEMENT  OF  A 
CONTINOUS  MEMBER  AT  DIFFERENT 

LEVELS  OF  Y    162 

H     CASE  4.   THE  AXIAL  STRESS  DISTRIBUTIONS 
OF  A  DOUBLE-DOUBLER  JOINT  AT 

DIFFERENT  LEVELS  OF  Y    171 

LIST  OF  REFERENCES    179 

BIOGRAPHICAL  SKETCH    183 


v 


LIST  OF  TABLES 

Table  Page 

3.1.  Elastic  material  properties  used  in 

analysis    33 

4.1.  Shear  stress  concentration  factor 

and  the  order  of  stress  singularity  x 

for  case  1  in  Chapter  5.3    82 

5.1.  Selected  material  properties  for 

singular  behavior  analysis    100 

5.2.  Different  values  Cy,         for  Case  1 
in  rx  direction  in  simple  tension 

loading  case    114 

5.3.  Different  values  Cx,  Cy,  C^  for  Case  2 
in  the  r2  direction  in  simple  tension 

loading  case    114 

5.4.  Different  values  C^  for  Case  1  in  the 
r3  direction  in  simple  tension 

loading  case    115 

5.5.  Different  values  Cx,  C^  for  Case  2 

in  the  rA  direction  in  simple  tension 

loading  case    115 


vi 


LIST  OF  FIGURES 


Figure  Page 
1.1  Various  types  of  joints    4 

2.1  The  single  lap  adhesive  bonded  joint    6 

2.2  Identification  of  the  various  parts  of 

a  single  lap  joint    6 

2.3  Systems  considered  by  Goland  and 

Reissner    10 

2.4  The  single  lap  joint  analyzed  by 

Pahoja    17 

3.1  Double-doubler  joint  configuration    31 

3.2  Double-doubler  joint  geometry    32 

3«3  Material  principal  coordinate  system    34 

3.4  Displacement  boundary  conditions  for 

simple  tension  and  pure  bending   38 

3«5  The  finite  element  mesh  coordinates 

for  the  high  stress  concentration 

regions    3g 

3.6  The  layout  of  finite  element  mesh  in 

simple  tension    41 

3.7  The  layout  of  finite  element  mesh  in 

pure  bending    42 

3-8             Case  1.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and  adhesive    44 

J. 9  The  different  levels  of  y  in  the  adhesive 

StriP    46 

MO  Case  1.  The  stress  distributions  of 

double-doubler  reinforcement  of  a  continous 
member  along  the  boundary  between 

the  outer  adherend  and  adhesive    48 

vii 


Case  1.  The  stress  distributions  of 
a  double-doubler  joint  along  the  boundary 
between  the  central  adherend  and 
adhesive   

Case  l.  The  stress  distributions  of 

a  double-doubler  joint  along  the  boundary 

between  the  outer  adherend  and 

adhesive   

Case  2 .  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and  adhesive   

Case  2.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  outer  adherend  and  adhesive   

Case  2.  The  stress  distributions  of 
a  double-doubler  joint  along  the  boundary 
between  the  central  adherend  and 
adhesive   

Case  2.  The  stress  distributions  of 

a  double-doubler  joint  along  the  boundary 

between  the  outer  adherend  and 

adhesive   

Case  3.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and  upper  adhesive  . 

Case  3.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  upper  adherend  and  adhesive  

Case  3.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and 
lower  adhesive   

Case  3.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  lower  outer  adherend  and  adhesive 


viii 


Case  3.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and 

upper  adhesive   

Case  3 .  The  stress  distributions  of 
a  double-doubler  joint  along  the  boundary 
between  the  upper  outer  adherend  and 
adhesive   

Case  3 .  The  stress  distributions  of 

a  double-doubler  joint  along  the  boundary 

between  the  central  adherend  and 

lower  adhesive   

Case  3.  The  stress  distributions  of 
a  double-doubler  joint  along  the  boundary 
between  the  lower  outer  adherend  and 
adhesive   

Case  4.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and 

upper  adhesive   

Case  4.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  upper  outer  adherend  and  adhesive  . . . 

Case  4.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and 
lower  adhesive   

Case  4.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  lower  outer  adherend  and  adhesive  . . . 

Case  4.  The  stress  distributions  of 
double-doubler  reinforcement  of  a 
continous  member  along  the  boundary  between 
the  central  adherend  and 
upper  adhesive   


ix 


3.3  0  Case  4.  The  stress  distributions  of 

a  double-doubler  joint  along  the  boundary 

between  the  upper  outer  adherend  and 

adhesive    75 

3.31  Case  4.  The  stress  distributions  of 

a  double-doubler  joint  along  the  boundary 
between  the  central  adherend  and 

lower  adhesive    76 

3.32  Case  4.  The  stress  distributions  of 

a  double-doubler  joint  along  the  boundary 

between  the  lower  outer  adherend  and 

adhesive    77 

4 • 1             Axial  lateral  normal  stress  distribution 
on  the  boundary  between  the  central 
adherend  and  adhesive  for  four  values 
of  E/Ea    83 

4.2  Axial  shear  stress  distribution 

on  the  boundary  between  the  central 

adherend  and  adhesive  for  four  values 

of  E/Ea    84 

4«3             Axial  lateral  normal  stress  distribution 
on  the  boundary  between  the  outer 
adherend  and  adhesive  for  four  values 
of  E/Ea    85 

4.4  Axial  shear  stress  distribution 

on  the  boundary  between  the  outer 

adherend  and  adhesive  for  four  values 

of  */*   86 

5.1  Two  edge-bonde  elastic  wedges  of 

different  materials  under  normal  and 
shear  loading   


5.2              Case  1.  Half  plane  bonded  to  a  quarter 
plane   


5.5  The  coarse  model  and  sub-model 
in  case  l   

5.6  The  coarse  model  and  sub-model 
in  case  2   


91 


99 


5.3  case  2.  Three  quarter  plane  bonded  to 

a  quarter  plane    104 

5.4  Specified  boundary  displacement  method  .. 


108 


112 


113 


x 


A.l             Case  1.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.06375    124 

A«2             Case  l.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  ■  0.0675    125 

A. 3             Case  1.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.08    126 

A.  4  Case  1.  Stress  distributions  of  double- 

doubler  reinforcement  of  a  continous 

member  at  the  level  of  y  =  0.0925    127 

B.  l  Case  1.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.06375    12g 

B.2             Case  1.  Stress  distributions  of 
a  double-doubler  joint  at  the 
level  of  y  =  0.0675   


B«3             Case  1.  Stress  distributions  of 
a  double-doubler  joint  at  the 
level  of  y  =  0.08   

B-4             Case  1.  Stress  distributions  of 
a  double-doubler  joint  at  the 
level  of  y  =  0.0925   


130 


131 


132 


C.l             Case  2.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.06375    134 

C.2              Case  2.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.0675    135 

C.3             case  2.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.08    13  6 

C.4              case  2.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.08    137 


xi 


Case  2.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.06375   

Case  2.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.0675   

Case  2.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.08   

Case  2.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.0925   

Case  3 .  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.095   

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.08   

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.07 

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.065   

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.065   

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.07 

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.08 

Case  3.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.095 

Case  3.  Stress  distributions  of 
a  double-doubler  joint  at  the 
level  of  y  =  -0.095   


xii 


Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.08   , 

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.07   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.065   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.065   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.07   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.08   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.095   

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.095 

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.08 

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.07 

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  -0.065 

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.065 

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.07 


xiii 


Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.08   

Case  4.  Stress  distributions  of  double- 
doubler  reinforcement  of  a  continous 
member  at  the  level  of  y  =  0.095   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.095   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.08   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.07   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  -0.065   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.065   

Case  3.  Stress  distributions  of 

a  double-doubler  joint  at  the 

level  of  y  =  0.07   

Case  3.  Stress  distributions  of 
a  double-doubler  joint  at  the 
level  of  y  =  0.08   

Case  3.  Stress  distributions  of 
a  double-doubler  joint  at  the 
level  of  y  =  0.095   


xiv 


Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfilment  of 
the  Requirements  for  the  Degree  of  Doctor  of  Philosophy 


STRESS  ANALYSIS  OF  A  DOUBLE-DOUBLER  JOINT 
IN  COMPOSITE  STRUCTURE 


BY 

MYUNG  KYUN  PARK 

May,  1990 

Chairman:  Lawrence  E.  Malvern 

Major  Department:  Aerospace  Engineering, 

Mechanics  and  Engineering  Science 

This  study  first  provides  a  linear  elastic  stress 
analysis  using  the  finite  element  method  for  a  double-doubler 
adhesive  joint.  The  joint  is  modelled  such  that  the  variation 
of  stresses  and  strains  through  the  adhesive  thickness  can  be 
determined.  The  linear  elastic  stress  analysis  is  focused  on 
three  areas:  first,  the  interfacial  stress  distribution 
between  the  adherend  and  adhesive;  second,  cohesive  stress 
distribution  in  the  adhesive;  and  third,  stress  distribution 
across  the  thickness  of  the  adhesive  layer  at  the  ends  of 
layer. 

The  second  step  is  to  carry  out  a  parametric  study  on 
the  double-doubler  joint.  since  the  stress  distribution  in 
the  adhesive  layer  depends  on  several  parameters  (geometrical 
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and  material) ,  an  understanding  of  the  influence  of  these 
various  parameters  will  lead  to  better  joint  design  for 
specific  requirements. 

The  third  step  is  to  devise  parameters  that  describe  the 
singular  behavior  near  the  bonding  edges  between  the  adherend 
and  the  adhesive  layer  by  applying  the  existing  elastic 
solution  for  stress  singularities  to  the  double-doubler  joint 
problem.  The  specified  boundary  displacement  method  is  then 
used  to  determine  the  values  of  parameters.  These  parameters 
can  be  used  as  a  means  of  failure  prediction  by  comparing  the 
calculated  parameters  with  known  critical  values  of  the 
parameters  for  the  material  and  comparison  of  stresses  very 
near  the  bonding  corners  where  stress  singularities  occur. 
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CHAPTER  1 
INTRODUCTION 


1 . 1  Introduction 
In  the  design  of  material  structures  components  must  be 
joined  in  such  a  manner  that  overall  the  structure  retains  its 
structural  integrity  while  performing  its  intended  function 
subjected  to  loads  and  environment.  Joining  metallic 
structures  is  a  well-developed  technology  involving  riveting, 
bolting,  welding,  glueing,  brazing,  soldering,  and 
combinations  thereof.  However,  for  polymer  matrix  fiber 
reinforced  composites  only  adhesive  bonding  and  mechanical 
fasteners  (bolts  and  rivets)  can  be  utilized. 

Recent    advances    in   composite    structures    and   the  new 
adhesive  bonding  techniques  based  upon  very  strong  epoxy  type 
adhesives  have  made  feasible  the  adhesive  joining,  stiffening 
and  repairing  of  structural   elements   subjected  to  extreme 
enviromental  and  loading  conditions.    This  type  of  joining  and 
stiffening    has    been    used    in    the    aerospace,  automobiles, 
electronics,      sporting     goods     and     appliance  industries. 
Consequently  the  importance  of  adhesive  bonding  in  technology 
has  been  recognized.      Various  types  of  adhesive  joints  are 
shown  in  Figure  1.1.     Among  the  various  types  of  joints,  the 
double-doubler  joint  is  selected  for  this  study. 
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Adhesive  bonding  provides  structures  having  lighter  weight, 
more  fatigue  resistance,  improved  aerodynamic  smoothness, 
better  distribution  of  imposed  stress,  better  thermal  and 
electrical  insulation  and  greater  crack  retardance  over 
mechanical  fasteners.  Therefore,  adhesive  bonding  is  very 
desirable  for  use  in  composite  material  structures. 

The  use  of  adhesives  is  growing  rapidly  but  designers 
lack  the  definitive  approach  that  gives  confidence  in 
routinely  incorporating  adhesive  bonding  into  design. 

1.2  Objective  and  Scope 
The  analysis  of  the  behavior  and  the  failure  of  an 
adhesive  joint  must  include  both  the  understanding  of  the 
mechanism  of  adhesion,  which  is  closely  related  to  the 
microstructure  between  the  adhesive  polymer  and  the  adherends, 
and  the  analysis  of  stresses  in  the  joint.  This  study  will 
be  focused  on  the  latter  part.  Failure  in  an  adhesive  joint 
can  be  either  "adhesive,"  which  occurs  at  the  interfaces 
between  the  adhesive  and  the  adherend,  or  "cohesive,"  which 
occurs  either  in  the  adhesive  or  in  the  adherends. 

Therefore  the  first  objective  of  this  thesis  is  to  do 
a  linear  elastic  stress  analysis  for  the  double-doubler  joint 
by  means  of  finite  element  method  (using  ANSYS) .  The  linear 
elastic  stress  analysis  is  focused  on  three  areas:  first,  the 
interfacial  stress  distribution  between  the  adherend  and  the 
adhesive;     second,     cohesive     stress     distribution     in  the 


adhesive;  and  finally,  the  stress  distributions  at  various 
levels  of  y  in  the  adhesive  layer  near  the  ends  of  layer. 

The  second  objective  is  to  carry  out  a  parametric  study 
on  the  double-doubler  joint.  Since  the  stress  in  the  adhesive 
layer  depends  on  several  parameters  (geometrical  and 
material) ,  an  understanding  of  the  influence  of  these  various 
parameters  will  lead  to  better  joint  design  for  specific 
requirements . 

The  third  objective  is  to  devise  parameters  that  describe 
the   singular   behavior   near   the   bonding   edges   between  the 
adherend    and    the    adhesive    layer    by    applying    an  existing 
elastic    solution    for    stress    singularities   to    the  double- 
doubler    adhesive    joint    problem.       The    specified  boundary 
displacement  method  (SBD)   is  used  to  devise  the  parameters. 
These  parameters  can  be  used  as  a  means  of  failure  prediction 
by  comparing  the  calculated  parameters  with  known  critical 
values  of  these  parameters  for  the  material  and  comparison  of 
stresses     very     near     the     bonding     corners     where  stress 
singularities  occur. 
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1.1  Various  types  of  joints 


CHAPTER  2 
REVIEW  OF  PREVIOUS  WORK 


2 . 1  Introduction 
A  considerable  amount  of  experimental  and  theoretical 
research  has  been  carried  out  on  adhesive  bonded  joints. 
Theoretical  and  experimental  stress  analyses  of  various  types 
of  adhesive-bonded  joints  have  been  made  by  many  authors.  Due 
to  the  relatively  simple  geometry  and  wide  application,  many 
theoretical  and  experimental  analyses  have  been  devoted  to  the 
single  lap  joint.    (See  Figure  2.1) 

In  this  review  the  important  investigations  of  the  single 
lap  joint  and  known  investigation  of  the  double  lap  joints  are 
included.  The  following  review  provides  a  history  and 
background  of  the  present  state  of  knowledge  in  adhesively 
bonded  joints.  First  the  investigations  of  the  single-lap 
joint  are  reviewed  here  because  the  earlist  analyses  and  the 
analyses  which  are  used  as  the  basis  for  most  subsequent 
analyses  were  made  on  the  single-lap  joint. 

The  terminology  used  in  identifying  various  parts  of  the 
joint  is  shown  for  a  single  lap  joint  in  Figure  2.2.  The 
"overlap"  is  the  length  over  which  the  adhesive  extends,  and 
the  "overhang"  is  the  part  of  the  adherend  outside  the 
overlap.     The  adhesive  layer  has  two  interfaces  and  two  free 


adherend  1 


adhesive 


Figure  2.1  The  single  lap  adhesive  bonded  joint 

l ead  i  ng  end 


^-overlap  -►^_overhang  _^ 


Figure  2.2  Identification  of  the  various  parts 
of  a  single  lap  joint 
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edges.  The  "leading  end"  of  an  adherend  is  where  the  overhang 
begins  after  the  overlap.  The  "trailing  end"  of  the  same 
adherend  is  at  the  opposite  end  of  the  overlap.  Corresponding 
parts  of  the  double-lap  joints  are  identified  with  the  same 
terminology.  A  typical  single-lap  joint  is  under  tensile 
loading. 

2.2  Analytical  Studies 
2.2.1  Volkersen ' s  Theory 

Classical  analysis  of  an  adhesive  bonded  joint  was  made 
by  Volkersen  [1].  He  derived  the  load  distribution  of  a 
multi-row  riveted  lap  joint  in  tension,  where  he  idealized 
the  rivets  by  replacing  them  by  a  continous  medium  of  given 
shear  flexibilty,  and  considered  the  case  of  an  adhesive  bond. 
Bending  of  the  adherends  and  the  associated  lateral  normal  and 
axial  normal  stresses  were  not  considered.  Therefore  the 
shear  stress  was  forced  to  be  constant  through  the  thickness 
of  adhesive.  By  setting  up  the  differential  equations  of 
linear  elasticity,  he  showed  that  shear  stress  at  any  point 
along  the  length  of  the  overlap  was  given  with  respect  to 
dimensionless  quantites: 

a)  the  overlap  ratio 

D  =  Gc2/E2tat2  21 

b)  the  adherend  stiffness 


The  shear  stress  distribution  in  the  adhesive  layer  was 
nonuniform  in  the  x-direction,  but  assumed  uniform  constant 
in  the  y-direction.  Peaks  of  the  shear  stress  occurred  at  the 
ends  of  the  overlap.  He  compared  the  maximum  shear  stress, 
rmax  at  the  ends  of  the  overlap  with  the  mean  stress,  rm  and 
found  the  stress  concentration  factor,  7,  ,  for  joints  with 
nonidentical  adherend. 


ii  = 


1/2      (W-l)cosh[  (Dwj1/2] 

=  (d/w)  1/2   .  ;  

sinh[  (DW)1/2] 


2.3 


where        D  =  Gc2/E2t2ta 

W  =  (Eiti+Eat^/Eiti 

rm  =  averaged  applied  shear  stress 

This  is  for  the  case  Ext,  >  E2t2.  He  showed  that  the  stress 
concentration  factor  increased  with  increasing  shear  modulus 
of  the  adhesive  and  with  decreasing  moduli  of  elasticity  and 
thickness  of  the  adherends. 

The   theory   developed   by   Volkersen   has   two  important 
shortcomings . 

a)  The  bending  deformation  was  neglected.  (There  should 
be  a  bending  moment  applied  to  the  joint  in  addition 
to  the  in-plane  tension,  since  the  two  applied  forces 
are  not  collinear.) 
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b)  The  assumption  was  made  that  the  shear  stress  does  not 
vary  across  the  thickness. 


2.2.2     Goland  and  Reissner  Theory 

Goland  and  Reissner  [2]  studied  the  stress  distribution 
for  adhesive  joints  with  identical  adherends  under  the  action 
of  tensile  forces  at  the  ends  of  the  adherends.     The  forces 
act  in  the  direction  parallel  to  the  bonding  line  before  any 
deformation  takes  place,  but  the  line  of  action  of  the  tensile 
force    R    is    assumed    to    pass    through   the   midpoint    of  the 
adhesive  layer  for  the  deformed  equilibrium  configuration. 
They  recognized  that  the  bending  of  the  adherends  beyond  the 
bonded    region    has    a    significant    effect    on    the  stress 
distributions    in   the    joint    itself.       They    approached  the 
problem  in  two  steps.     First  they  determined  the  loads  at  the 
edges  of  the  joint  which  were  transmitted  from  the  loads  at 
the  ends  of  the  adherends.     In  the  first  step  of  the  problem, 
they  treated  the  adherends  as  cylindrically  bent  plates,  and 
found   that    forces    transmitted   to   the    edges    of   the  joint 
consisted  of  a  tensile  force  P,  a  small  shear  force  V,  and 
bending  moment  M,  as  shown  in  Figure  2.3  in  which 


a 


P  =  R 

>n  n  -  ,,  2\  t>/i?  4.  -i  1/2 

2.4 


V  =  KP[3(1  -  i^JP/EitJ1'2 


M  =  KPtx/2 


Figure  2.3  Systems  considered  by  Goland  and  Reissner 
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where  K  is  the  bending  moment  factor  which  is  the  ratio  of  the 
magnitude  of  the  edge  moment  M  and  the  value  of  this  moment 
in  a  nondeformable  system: 


K  1   2.5 

1+  (2/2)  tanh[  (c/itj  J&{l-v\)  P/Eiti  ] 


The  value  K=l  indicates  that  the  adherends  exhibit  no 
deformation  due  to  bending,  and  K  <  1  indicates  that  a  certain 
amount  of  bending  deformation  takes  place. 

In  the  second  step  of  the  problem,  they  assumed  that  the 
normal  stress  in  the  adhesive  parallel  to  the  layer  can  be 
neglected,  and  that  the  lateral  normal  stress  and  shear  stress 
do  not  vary  across  the  thickness  of  the  adhesive  layer.  They 
used  a  plane  strain  analysis  and  solved  the  problem  for  two 
limiting  cases: 

1)  The  case  where  the  adhesive  layer  is  so  thin  and  stiff 
that  its  deformation  may  be  neglected.  This  case  requires 
that  the  condition  ta/Ea  «  tl/El.  The  results  show  that, 
while  the  lateral  normal  stress  is  high  at  the  edge  of  the 
joint,  the  shear  stress  is  zero  there.  The  shear  stress, 
however,  rises  rapidly  to  a  maximum  value  in  close  proximity 
to  the  edge. 

2)    The   case   in  which  the   adhesive   layer   is   soft  and 
flexible    and   the    joint    flexibility    is   mainly   due   to  the 
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deformation  of  the  adhesive  layer.  This  case  requires  the 
condition  t./Ea  »  t^E^  The  shear  stress  r  and  lateral 
normal  stress  a  are  expressed  as  follows. 


T  = 


P         r     Be  cosh(fix/t,l 

4c     I  TET  (1+3K)  coBhcflcAti)  +  3(1_K)  ]  2-6 


a  =(1/A)  (P/tJ  (2^0)2   [  (R2A2  K/2  + 

AK'coshA  cosA) cosh(2Ax/c)   cos(2Ax/c)  2.7 
+  (RiA2K/2  +  AK'sinhA  sinA ) sinh (2Ax/c)  sin(2Ax/c)] 
where 

B  =     (SGtJ  Ext)1/2 

A  =  fe/2ta)  (SEti/Eit)1'* 

R:  =  coshA  sinA  +  sinhA  cosa 

R2  =  sinhAcosA  -  coshAsinlA 

a  =  (1/2) (sinh2A  +  sin2A) 

K«=  Vc/2Pt! 

The  maximum  value  of  r  and  a  are  found  at  the  edges  of  the 
joint. 

fP/tJ  (    t,  )      -     T   t  -2t[   (1+3K)    coth  -ff-  +  3(1-K)    ]  2.8 


To/Ti     <  It   >2   ■  -Si"h?A  ~  S.in2A    +  \jri   CQSh2A  +  soggA 

(P/tJ      2t,  2         smh2A  +  sin2A    +         sinh2A  +  sin2A 


2.9 
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Compared  with  the  factor  D  used  by  Volkersen,  it  can  be  seen 
that  &c/2tx  =  (2D)1/2.  For  large  values  of  D,  coth(2D)1/2 
approaches  the  value  of  1,  and  therefore  the  shear  stress 
concentration  factor,  rll,  can  be  written,  for  large  D,  as 

1l  =  rmax/(p/c)    «   (1/4)  (1+3K)  (2D)1/2  +  3(l-K)/4  2.10 
where  D  =  Gc2/E2t2ta 

The  limitations  of  Goland  and  Reissner  theory  include 

a)  It  is  valid  only  for  identical  adherends. 

b)  it  is  valid  only  for  two  special  cases;  the  first  case 
is  for  VE,  «  ti/Ei,  while  the  second  case  reguires 
that  tJEa  »  taEx. 

c)  The  stress  across  the  adhesive  layer  thickness 

is  assumed  to  be  constant,  which  results  in  inaccurate 
predictions  in  the  region  adjacent  to  the  edges  of  the 
overlap. 


2.2.3  Plantema's  Analysis 

Plantema  [3]  combined  the  Goland-Reissner  theory  with 
Volkersen -s  theory  in  the  analysis  of  a  single-lap  adhesive 
bonded  joint.  The  factor  K  from  the  Goland-Reissner  theory 
was  employed  to  calculate  the  differential  strain  of  the 
members  due  to  bending,  while  Volkersen -s  theory  was  used  to 
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calculate  the  differential  strain  of  the  members  due  to  shear, 
and  the  stress  at  the  edges  of  the  overlap. 

The  shear  stress  concentration  factor  was  given  as 

91  =   [D(l+3K)/2]1/2coth[D(l+3K)/2]1/2  2.11 

He  did  not  make  an  attempt  to  calculate  the  lateral  normal 
stress. 

2.2.4  Kelsey  and  Benson  Theory 

Kelsey  and  Benson  [4]  employed  the  complementary  energy 
method  to  determine  the  shear  and  lateral  normal  stresses  in 
a  single  lap  joint.  A  linear  variation  of  the  lateral  normal 
stress  was  allowed  across  the  the  adhesive  layer,  while  the 
shear  stress  was  assumed  constant  across  the  adhesive 
thickness.  The  result  was  a  differential  equation  of  higher 
order  than  that  obtained  by  Volkersen,  which  enabled  the 
boundary  condition  of  zero  shear  stress  at  the  ends  of  the 
overlap  to  be  satisfied.  The  shear  and  lateral  normal  stress 
obtained  are  as  follows: 
Let  A  =  pm^c 

D  =  (%coshmlC/2)  (sinhm2c/2)-(m2sinhm1c/2)  (coshm2c/2) 


r  -  (A/4D)  [(cosh^x)  (coshm2c/2)-(coshmlC/2)  (coshm2x)  ] 

2.12 


(At/8D)  [  (niiSinhmiX)  (coshm2c/2) 
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-(m2coshm1c/2)  (sinhm2c/2)  ] 

2.13 


The  shear  stress  concentration  factor  r)  =  r/(p/c)  is  given  by 

v  =    m^c^coshro,  (xc/2) )  (coshm-c/2)  (coshm,c/2)  fcoshm,  (xc/2) 
4  m1(coshm1c/2)  (sinhm2c/2) -m2(sinhm2c/2)  (coshm2c/2) 

2.14 


where 


m2    =  Gt'» 


'a 


m3  -  (  i  -        2G2fcB  )  V2 

Kelsey  and  Benson  imposed  the  zero  shear  stress  boundary 
condition  which  causes  a  sharp  drop  of  shear  stress  from  its 
peak  value  at  positions  very  near  the  ends  to  zero  at  ends. 
They  reasoned  that  this  is  necessary  for  equilibrium  because 
the  adhesive/air  boundary  can  not  sustain  a  shear  stress. 

2.2.5  Pahoia's  Analysis 

Pahoja  [5]  made  the  analysis  of  the  joint  under  tensile, 
shear,  and  bending  loads  as  shown  in  Figure  2.4.  The  joint 
was  subjected  to  a  general  loading  condition.  He  employed  the 
theory  of  minimum  potential  energy  to  study  the  axial  normal, 
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lateral  normal  and  shear  stresses  in  the  adhesive  layer.  He 
defined    the    displacements    in    the    adhesive    in    terms  of 
polynominals  with  m  terms   for  the  displacements   in  the  x 
direction  and  n  terms  for  the  displacement  in  the  y  direction. 
The  total  energy  was  expressed  in  terms  of  the  displacement 
components  in  series  form.    By  minimizing  the  total  potential 
energy  of  the  system,   a  system  of  governing  linear  ordinary 
differential  equations  and  natural  boundary  conditions  were 
obtained.     The  equations  were  reduced  to  algebraic  and  first 
order  differential   equations  by  substituting  new  variables 
for  the  displacements  and  their  derivatives.    By  solving  these 
equations   numerically,    the   stress   distribution   and  stress 
concentrations  were   found.      The   stress   concentrations  for 
identical  adherends  have  the  follwing  forms: 

Shear  stress 

('maxAav)     =    1    +     ( 6/7 )[( 19/18 )     -    Vf]  (Sf    -  1) 

for  Vf  <  l  and  Sf  >  4 

2.15 

('maxA.v)     =    "1    +     (18/19)  (Vf    -    17/18)  (Sf    -  1) 

for  Vf  >  1  and  Sf  >  5 


Lateral  normal  stress 

K^axAav   =    (t2/c)(1.4    -    1.675    Vf)  Nf2 


2.16 


17 


Figure  2.4  The  single  lap  joint 
analyzed  by  Pahoja 


18 


Axial  normal  stress 
(OnmxAav  =    (ti/c)(0.85  -  0.95V£)    Nf2  2.17 

where 

V   c 


Vf  = 


T         [t.  +  Q,  t/  ft  +  M/fTyj 


Sf  =   (c/tj  (Et^Et.)1'2 
Nf  =   (C/tJ  (Et^Eit,)1'* 

Qi   -  Ei/fl-i/^) 

*  =  Qt.2/6  +  Qjtjt,  +  Qiti2 


Q  =  E/(l-„2) 


Pahoja  concluded  that  the  single  lap  joint  with  identical 
adherends  was  the  best  design.  This  joint  was  most  efficient 
when  loaded  such  that  the  shear  stress  distribution  was 
uniform,  i.e.,  vf  =  l.  He  compared  his  theoretical  results 
with  the  results  of  photoelastic  analysis,  and  these  are  in 
good  agreement.  The  limitation  of  this  analysis  is  that  the 
free  boundary  condition  at  the  ends  of  the  adhesive  were 
violated  just  as  in  Volkersen's  and  Goland  and  Reissner's 
analyses. 
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2.2.6  Renton  and  Vinson's  Analysis 

Renton  and  Vinson  [6,7]  made  the  most  complete  analysis 
of  a  single  lap  joint.  They  clearly  identify  the  extent  of 
the  influence  of  the  various  parameters  on  the  stress 
distribution  in  the  adhesive  and  adherends.  The  analysis  was 
very  general  and  applicable  to  isotropic  and  anisotropic 
adherends  as  well  as  to  dissimilar  and  similar  adherends.  The 
method  of  analysis  that  they  developed  involve  the  following 
assumptions  and/or  limitations: 

1)  The  laminated  adherends  are  symmetric  about  their  own 
midsurfaces  (i.e.,  no  bending  stretching  coupling). 

2)  Each  ply  or  lamina  in  each  adherend  is  orthotropic 

3)  The  effective  elastic  mechanical  properties  of  the 
adhesive  are  accounted  for. 

4)  Both  shear  and  lateral  normal  stresses  are  accounted 
for  in  the  adhesive;  they  vary  in  the  load  direction, 
but  do  not  vary  in  the  thickness  direction. 

5)  The  adhesive  thickness  is  much  smaller  than  the 
adherend  thickness,  i.e.,  ta  «  t,  and  ta  «  t2 

6)  Shear  deformation  and  lateral  normal  strains  are 
accounted  for  in  each  adherend. 

Composite  material  laminated  plate  theory  was  used  as  the 
building  block  or  starting  point  for  developing  the  method  of 
analysis.  The  lap  joint  structure  was  divided  into  five 
portions  for  analysis:  the  adhesive,  the  adherends  on  either 
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side  of  the  adhesive,  and  the  adherend  portions  from  the  end 
of  lap  joint  to  the  load  application.  Twenty-six  boundary 
conditions  are  specified  and  satisfied. 

The  analytical  solutions  obtained  for  stresses  and 
deformations  throughout  the  structure  were  programmed  on 
computer  for  easy  calculation.  After  developing  a  valid 
closed  form  solution  for  the  single  lap  joint  problem,  they 
ascertained  the  influence  of  certain  important  parameters  on 
the  stress  distributions  in  the  adhesive. 

The  influence  of  the  joint  parameters  on  the  stress 
distributions  in  the  adhesive  are  summarized  below: 

1)  The  greater  the  difference  in  the  moduli  of  the 
adherends,  the  higher  is  the  magnitude  of  the  maximum 
shear  stress. 

2)  For  dissimilar  adherends,  the  maxima  of  both  the 
shear  and  the  lateral  normal  stresses  in  the  adhesive 
occured  at  the  leading  end  of  the  less  rigid  adherend. 

3)  The  greater  the  ratio  of  the  moduli  of  the  adherend 
and  the  adhesive,  the  more  uniform  the  shear  stress 
distribution. 

4)  The  larger  the  ratio  of  the  moduli  of  the  adherends 
and  the  adhesive,  the  smaller  is  the  magnitude  of  the 
maximum  lateral  normal  stress. 
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In  the  design  of  single  lap  adhesive  bonded  joints,  they 
recommended  that 

1)  The  adherends  of  joints  should  be  identical. 

2)  The  modulus  of  the  adherends  should  be  high  compared 
to  that  of  adhesive. 

3)  The  length  of  overlap  should  be  about  ten  times  the 
minimum  thickness  of  the  adherend. 

4)  The  joint's  intended  function  should  influence 
adhesive  selection:  in  static  bodies,  adhesive 
ultimate  strength  should  be  high;  in  fatigue,  adhesive 
fracture  toughness  should  be  considered. 

2.2.7  Hart-Smith's  Analysis 

The  analysis  of  a  double  lap  joint  by  Hart-Smith  [8]  was 
designed  to  calculate  the  joint  efficiency  of  a  double  lap 
joint  in  terms  of  the  strength  of  the  joint. 

The  adhesive  was  characterized  as  elastic-plastic  in 
shear.  In  the  analysis,  the  shear  stress  distribution  was 
constant  across  the  adhesive  thickness. 

In  his  analysis,  the  deformation  of  the  adhesive  was 
greatest  at  the  ends  of  the  overlap,  and  consequently  the 
behavior  of  the  adhesive  was  characterized  as  plastic  at  the 
ends  of  the  lap. 

In  the  elastic  analysis,  the  shear  stress  in  the  adhesive 
was  expressed  in  terms  of  the  displacements  in  the  outer  and 
center  adherends.      The  differential  equation  governing  the 
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shear  stress  distribution  was  determined  from  the  stress 
strain  relationship  and  the  force-equilibrium  equation. 

In  the  elastic-plastic  analysis,  the  plastic  shear  stress 
and  the  distance  over  which  the  adhesive  was  plastic  were 
determined  from  the  boundary  condition  for  the  plastic  region, 
the  equation  ensuring  the  continuity  of  the  elastic  and 
plastic  shear  strains,  and  the  gross  horizontal  equilibrium 
equation  of  the  joint.  The  expression  for  the  plastic  shear 
stress  was  then  rearranged  to  express  the  applied  load  in 
terms  of  the  shear  strain  energy  per  unit  area  of  the  adhesive 
layer. 

The  lateral  normal  stress  in  the  adhesive  was  considered 
to  be  independent  of  the  shear  stresses.  The  analysis  for  the 
lateral  normal  stress  was  applicable  only  when  the  adhesive 
was  deformed  beyond  the  elastic  range  in  joints  with  long 
overlaps.  The  lateral  displacement  of  the  outer  adherends 
was  related  to  the  moment  at  the  ends  of  the  overlap  by  the 
theory  of  bending  of  plates. 

The  coefficients  of  the  normal  displacements,  which  were 
expressed  in  terms  of  trigonometric  and  hyperbolic  functions, 
were  determined  from  the  boundary  conditions  in  terms  of  the 
maximum  shear  stress.  The  lateral  normal  stress  was  then 
given  in  terms  of  an  "effective"  tension  modulus  of  the 
adhesive,  whose  value  is  difficult  to  establish  theoretically. 
The  conclusions  drawn  by  Hart-Smith  on  the  stress 
distributions  are 
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1)  The  highest  shear  stress  in  the  adhesive  was  at  the 
end  of  overlap  near  the  leading  end  of  the  less  stiff 
adherend . 

2)  For  long  overlaps,  the  magnitude  of  the  maximum  shear 
stress  did  not  vary  with  length;  only  the  stress  in 
the  central  portion  of  the  joint  was  lowered  by 
increasing  length. 

3)  Lateral  normal  stresses  were  reduced  by  using  a  thin 
outer  adherend,  or  an  adherend  with  a  low  extensional 
modulus,  or  a  thick  adhesive  layer.  Thus,  the  lateral 
normal  stresses  can  be  reduced  at  the  expense  of  the 
shear  stresses. 

4)  The  characteristic  extent  of  the  end  zones  was 
independent  of  the  length  of  overlap. 

2.3  Finite  Figment  Sfcyfllfeg 
2.3.1  Introduction 

The  use  of  finite  element  methods  of  analysis  means,  in 
principle,  that  the  simplifying  assumptions  made  in  the 
classical  analytical  methods  need  not  be  made.  The  technique 
does  allow  the  solution  of  problems  that  are  totally 
intractable  by  classical  methods.  Basically,  one  can  choose 
between  a  large  number  of  primitive  elements,  i.e.,  elements 
in  which  strain  is  assumed  to  be  constant,  or  a  smaller  number 
of  more  sophisticated  elements,  ones  with  linear  or  quadratic 
strain  variation.    Whatever  idealization  is  adopted,  a  higher 
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density  (finer  mesh)  of  elements  is  necessary  in  regions  of 
rapidly  varying  stress. 

Finite  element  analysis  of  an  adhesive-bonded  joint  was 
first  performed  in  1969.  Considerable  work  has  been  done  with 
finite  element  analysis  of  adhesive-bonded  joints  of  fiber 
reinforced  composite  materials. 

2.3.2  Ahluwalia's  Analysis 

The  first  finite  element  analysis  of  an  adhesive  bonded 
joint  was  made  by  Ahluwalia   [9].      He  analyzed  a  double-lap 
joint.     The  finite  element  solution  for  the  double  lap  joint 
was  for  plane  strain  with  120  nodal  points  and  188  triangular 
elements.     The  adhesive  was  divided  into  ten  equal  elements 
along   its   length   and  was   made   of   one   element   across  the 
thickness.     Consequently,  the  high  magnitude  of  the  stresses 
at   and   near   the   ends   of   the   overlap   were   not  correctly 
determined,  and  the  stresses  did  not  vary  across  the  thickness 
of  the  adhesive.     Smaller  elements  were  not  defined  at  the 
ends  of  the  overlap  where  shear  and  lateral  normal  stress 
concentrations  exist.     The  shear  and  lateral  normal  stress 
distributions  were  nonuniform  with  peaks  at  the  ends  of  the 
overlap.      Large  sizes  of  elements  in  areas  of  very  sharply 
rising    peak    stresses    result    in    lower    magnitude    of  the 
calculated  stresses. 
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2.3.3  Wooley  and  Carver's  Analysis 

Wooley  and  Carver  [10]  analyzed  a  single  lap  adhesive- 
bonded  joint  with  the  finite  element  method.  The  finite 
element  solution  for  the  single  lap  adhesive-bonded  joint  was 
for  plane  stress  problem  with  497  nodal  points  and  394 
elements. 

The  solution  was  based  on  the  computer  program  "ARBAS" 
by  E.L.  Wilson  of  the  University  of  California  at  Berkeley. 
The  locating  of  nodal  points  was  accomplished  by  dividing  the 
configuration  into  three  parts,  the  section  of  the  lap  joints, 
and  the  sections  to  each  side  of  the  lap.  in  the  program, 
quadrilateral  elements,  composed  of  four  triangular  elements, 
were  used  to  describe  the  joint.  The  displacements  of  the 
triangular  elements  were  described  by  linear  functions  which 
resulted  in  a  constant  strain  and,  therefore,  constant  stress 
over  each  element. 

The  adhesive  layer  was  divided  into  two  equal  elements 
through  its  thickness  and  twenty  equal  elements  along  its 
length.  The  maximum  aspect  ratio  of  a  quadrilateral  element 
in  the  adhesive  was  100:1.  The  single-lap  adhesive-bonded 
joint  was  made  of  identical  adherends  with  an  elastic 
adhesive.  The  Poisson's  ratios  of  the  adhesive  and  adherends 
were  0.3  for  all  the  cases.  The  range  of  the  ratio  of  the 
elastic  moduli  of  the  adherend  and  adhesive,  E/Ea  was  from  0.1 
to  1000.  The  range  of  the  ratio  of  the  thickness  of  the 
adhesive  and  adherends,  tjt  was  from  0.02  to  0.1.    The  length 
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of  overlap  was  taken  as  5,  10  and  20  times  the  adherend 
thickness. 

The  stress  distributions  at  the  adhesive-adherend 
interface  were  not  given.  The  stress  distributions  at  the  mid 
plane  of  the  adhesive  layer  were  nonuniform  with  peaks  near 
the  ends  of  the  overlap  for  both  shear  and  lateral  normal 
stresses.  They  compared  the  stress  concentrations  obtained 
by  each  method  at  three  values  of  the  geometric  and  material 
parameters,  with  two  elements  across  the  adhesive  thickness, 
the  variation  of  the  stresses  across  the  adhesive  thickness 
can  not  be  adequately  represented. 

2 . 4  Experimental  Studies 
2.4.1  Mclaren  and  Maclnnes's  Investigation 

Mclaren    and    Maclnnes    [11]    conducted  photoelasticity 
experiments    on    single    lap    joints.       The    object    of  their 
investigation   was   to   verify   the   mathematical    analysis  of 
Goland  and  Reissner.    Two  series  of  tests  were  performed,  in 
the  first  series,   the  complete  joint  was  made  of  Araldite. 
In  the  second  series  of  tests,    the  adherends  were  made  of 
aluminum  and  the  adhesive  of  Araldite.      The  ratio  of  the 
modulus  of  the  adhesive  to  that  of  the  adherend  was  1/20  in 
the  second  series  of  tests.     The  results  were  presented  in 
terms    of    the    fringe    order,    and    no    attempt    was    made  to 
determine   the  magnitude   of  the   stresses.      The  difference 
between  these  two  model  systems  was  in  the  transfer  of  load 
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from  one  adherend  to  the  other.  The  adhesive  in  the  second 
series  was  in  longitudinal  shear  across  its  entire  thickness; 
thus,  the  transfer  of  load  was  primarily  by  shear.  Their 
investigation  supports  the  theory  of  Goland  and  Reissner  in 
that  the  qualitative  aspects  of  the  theory  are  well 
corroborated . 


2.4.2  Tuzi  and  Shimada's  Investigation 

Tuzi  and  Shimada  [12]  conducted  a  photoelastic 
investigation  of  a  single  lap  joint.  Their  investigation 
concentrated  mainly  on  stress  concentration  in  the  joint. 
Single  lap  joint  models  having  different  adhesive  shapes  at 
the  ends  of  overlap;  rectangular  type,  with  a  fillet  and  with 
convex  and  concave  semi-circular  ends  were  used  in  the 
photoelastic  investigation. 

The  adhesive  layer  of  the  joint  was  made  of  epoxy  rubber 
and  the  adherends  were  represented  by  epoxy  resin  and 
duralumin.  In  the  first  experiment,  the  adherends  were  made 
of  epoxy  resin  KT-102  and  the  adhesive  layer  was  made  of  an 
epoxy  rubber.  The  adhesive  was  bonded  to  the  adherends  by 
Araldite  D.  in  the  second  experiment,  the  adherends  were  made 
of  duralumin,  the  adhesive  layer  was  made  of  epoxy  resin  KT- 
102,  and  the  adhesive  was  bonded  to  the  adherends  by  Araldite 
D.  Their  experimental  results  indicated  that  the  joint  with 
concave  semicircular  ends  was  the  best,  the  one  with 
rectangular  ends  was  second  best,  and  the  one  with  convex 


28 

semicircular  ends  was  the  worst  as  regards  the  stress 
concentration  in  the  joint.  They  also  studied  the  stress 
distributions  in  scarf  and  butt  joints  photoelastically . 
Their  photoelastic  models  were  made  of  epoxy  resin  and  epoxy 
rubber  and  also  epoxy  resin  and  metal.  Tension  and  bending 
tests  were  carried  out  on  bars  bonded  at  various  angles.  For 
the  joints  with  a  soft  adhesive  layer  bonded  to  hard 
adherends,  the  stress  concentrations  were  maximum  at  the  ends 
of  the  bonded  boundary  with  an  obtuse  angle.  For  the  joints 
with  a  hard  adhesive  layer  bonded  to  soft  adherends,  the 
stress  concentrations  were  found  at  the  ends  of  the  bonded 
boundary  with  an  acute  angle. 


CHAPTER  3 
LINEAR  ELASTIC  STRESS  ANALYSIS 


3.1  Introduction 

The  joint  under  consideration  is  a  double-doubler  joint 
(Figure  3.1).  This  joint  is  subjected  to  two  different 
loading  cases  (constant  tensile  loading  and  pure  bending) . 
In  this  chapter  the  three  materials  (central  adherend, 
adhesive,  outer  adherend)  are  assumed  to  behave  in  a  linear 
elastic  manner  under  the  applied  loads.  The  elastic 
properties,  Young's  modulus  and  the  Poisson's  ratio  are  Ec  and 
j/c  for  central  adherend;  Ea  and  »/a  for  the  adhesive  layer;  E0 
and  i/Q  for  the  outer  adherend.  The  material  properties  are 
extended  to  orthotropic  material  properties  for  the 
application  to  composite  adherends. 

The  overall  problem  is  to  find  the  stress  distributions 
in  the  double-doubler  adhesive  bonded  joint  under  the  applied 
load.  The  distributions  along  the  axial  direction  of  the 
stresses  on  the  interface  between  the  central  adherend  and 
adhesive,  and  of  the  cohesive  stresses  at  various  levels  in 
the  adhesive  are  the  particular  interest  in  this  study.  The 
analysis  can  be  simplified  to  either  a  plane  stress  or  a  plane 
strain  problem.  For  most  practical  purposes  the  width  of  an 
adhesive   joint   is   large  compared  to  the  thickness   of  the 

29 


30 

adherend    and    the    adhesive.       Therefore    the    plane  strain 
analysis   option   is  used  throughout.      First  of   all  linear 
elastic  stress  analysis  for  double-doubler  reinforcement  of 
a   continuous   member    (Figure   3.1a)    was   carried   out.  The 
results,  are  compared  with  the  previous  experimental  results 
and  they  show  good  agreement.     After  that  the  linear  elastic 
stress  analysis  for  a  double-doubler  joint  (Figure  3.1b)  was 
investigated.     The  finite  element  method  of  stress  analysis 
is    applied    to    this    plane    problem     involving  different 
materials.    The  reason  why  the  finite  element  method  is  chosen 
over  any  other  method  is  that  geometrically  complex  bodies  of 
different  materials  can  be  easily  represented  and  arbitrary 
distributions  of  displacement  and  stress  boundary  conditions 
can  be  specified. 

For  this  linear  elastic  stress  analysis  the  finite 
element  program  ANSYS  is  employed  in  the  finite  element 
studies  of  the  double-doubler  joint. 

3.2  Joint  Description 
The  geometry  and  dimension  of  the  joint  are  shown  in 
Figure  3.2.    The  basic  material  properties  used  in  this  joint 
are  given  in  Table  3.1. 
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a}  the  doub I e-doub I er  reinforcement 
of  a  continous  member 


b}  the  doub I e-doub ler  joint 


Figure  3.1  Double-doubler  joint  configuration 
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Figure  3.2  Double-doubler  joint  geometry 
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Table  3.1.  Elastic  Material  Properties  Used  in  Analysis 


Mateial 

Modulus , E 

Poisson's  ratio 

Shear  Modulus, G 

aluminum 

10.6E6  psi 

0.33 

epoxy 

3.6E6  psi 

0.39 

graphite 
epoxy 

22.8E6"  psi 
1.58E6b  psi 

0.3C 

0.81E6d  psi 

a:  Longitudinal  modulus  of  unidirectional  composite  i.e, 

fiber  direction 
b:  Transverse  modulus  of  composite 

c:  Major  Poisson's  ratio  of  unidirectional  composite 
d:  Inplane  shear  modulus 


3 . 3  Determining  Material  Properties 
The  strain-stress  relations  for  anisotropic  materials  can 
be  expressed  as 

€i  =  su  °i  ifj  =  1,2,  ,6  3.1 

where  £i  are  the  strain  components,  Sij  are  the  elements  of 
the  compliance  matrix,  and  o3  are  the  stress  components. 
Equation  3.2  gives  the  form  of  the  relations  for  a 
unidirectional  fiber  reinforced  lamina  with  axes  xlf  x2,x3 
aligned  with  the  principal  material  directions  as  shown  in 
Figure  3.3.  (xx  corresponds  to  fiber  direction,  x2  is  the  in- 
plane  coordinate  perpendicular  to  xlt  and  x3  is  the  tranverse 
coordinate  which  is  perpendicular  to  the  Xj-x2  plane.) 
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°\ 
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^23 

s33 
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0 

0 

°3 

0 
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0 

0 

°k 

0 

0 
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0 
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0 

*6 
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0 

0 

0 

^66 

3.2 


Figure  3.3  Material  principal  coordinate  system 
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The  compliances  can  be  expressed  in  terms  of  engineering 
constants  as  follows: 
Sn  =1/En 


3.3 


Si2  =  -^12/Eu  S13  =  -*/13/Eu 

S22  =  1/E22  S23  =  -u23/E22  S33  =  I/E33 

SA4  =  1/G12  S55  =  1/G23  S66  =  1/G13 

where 

Eiif  E22,  E33  are  Young's  moduli  in  1,2  and  3  directions 
respectively. 

i/ij  is  Poisson's  ratio  determined  from  contraction 
in  the  Xj  direction  during  a  tensile  test  in  the  xA 
direction 

Gl3  is  shear  moduli  in  the  x^Xj  planes 


We  know  that  the  following  relationships  exist  between 
the  Young's  moduli  and  the  Poisson's  ratios  due  to  symmetry 
of  the  compliance 

&U*21  =  E22»/12/  E22i/32  =  E33i,23,  E33l/13  =  Eul/31  3.4 
Thus  for  an  orthotropic  material,  constitutive  relations  are 


«2 


«  5 
«6 


1/En 

-^12/Eu 

-♦'is/En 
0 
0 
0 


1/E22 

_1/23/E22 

0 
0 
0 


-"la/Eu 

_,/23/E22 

VE33 
0 
0 
0 


0  0  0 

0  0  0 

0  0  0 

1/G12  0  0 
0  1/G23  0 
0  0  1/G13 


3.5 
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For  a  transversely  isotropic  material  in  which  the  x2-x3 
plane  is  the  plane  of  isotropy,  we  find  E22  =  E33,  G13  =  G12,  vn 
=  i/13  and  G23  =  E22/2  (1+j/23)  .  If  we  invert  Equation  3.5,  for  a 
material  having  tranverse   isotropy  relative  to  the  x,  -  x 

2  3 

plane,  we  obtain  the  elastic  stiffness,  C^,  in  terms  of 
engineering  constants  with  the  following  results: 

Cu  =  (1  -»223)  Ejj/V 

ci3  -  C12  =»/12(l  +  i/23)  E22/V 

C23  =   (f23  +  "2i2E22/Eu)  E22/V 

C33  =  C22  =  (1  -  v\2E22/En)  E22/V  3>6 
c*a  =  G23  =  E22/2(l  +  »/23) 

C55    =    C66    =  G12 

where  V  -  [(1  +  „„)  (1  -  „23  -  2„212E22/En)  ] 
For  the  plane  strain  case,Eu,   E22,   E33,    ,12,    „„,  Cuf  ^ 

G23  would  be  used  as  input  data  for  ANSYS .  If  G12,G13/G23  are 
not  specified,  they  are  computed  as: 

G12  =  E1E2/(E1  +  E2  +  2«/21Bi) 

G23  =  E2E3/(E2  +  E3  +  2i/32E2)  3<? 

G13  =  EiEa/fEi  +  E3  +  2*13E1) 

3.4  Finite  Element  Modgl  Represenl-aHnn 
Only  one  quarter  of  the  model  is  selected  for  analysis 
of  the  simple  tension  loading  case,  because  it  has  double 
symmetry  in  geometry,  materials,  boundary,  and  loading 
conditions.  One  half  of  the  model  is  chosen  for  the  pure 
bending  loading  case. 
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The     ANSYS     two-dimensional     isoparametric     element  which 
accomodates   orthotropic  material   properties   is   used.  The 
element   is   defined   by   four  nodal   points,    each   having  two 
degree  of  freedom.      The  plane  strain  option  is  used.  The 
displacement  boundary  conditions  imposed  on  simple  tension  and 
pure  bending  case  are  shown  in  Figure  3.4.     In  simple  tension 
cases   analyzed   in  this   study,    the  double-doubler  joint  is 
subjected  to  50  lb  at  the  outer  end  of  the  central  adherend. 
In    pure    bending    case,     a    linear    varying    force  boundary 
condition  equivalent  to  a  couple  is  used  at  the  outer  end  of 
the   central   adherend.       For   the    input   of   the  orthotropic 
material  properties  the  classical  thin  plate  lamination  theory 
is  used  to  calculate  the  effective  material  properties  for  the 
laminate.    For  unidirectional  composite  lamina,  properties  in 
the  thickness  direction  are  calculated  assuming  that  E22  =  e33, 
"12  =  "is  and  j/23  =  matrix  Poisson's  ratio. 

The  material  principal  coordinate  system  (1,2,3)  refers  to 
the  fiber  and  two  transverse  matrix  directions,  respectively. 
(See  Figure  3.3)  The  computer  program  can  be  used  to 
calculate  effective  material  properties  for  composite 
laminates.  The  mesh  is  so  generated  that  around  the  bonding 
corner  and  on  the  boundary  between  two  different  materials  the 
grids  are  much  finer  than  those  of  other  locations  for  better 
results  in  these  regions  of  high  stress  gradients.  The  finite 
element  dimensions  which  are  used  in  the  analysis  for  the  high 
stress  concentration  regions  are  shown  in  Figure  3.5. 
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Figure  3.4  Displacement  boundary  conditions  for 
simple  tension  and  pure  bending 
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Figure  3.5  The  finite  element  mesh  coordinates  for  the  high 
stress  concentration  regions 
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A  total  of  1591  nodes  and  1456  elements  for  the  simple  tension 
case,  and  2819  nodes  and  2672  elements  for  the  pure  bending 
case  are  used.  (See  Figure  3.6  and  Figure  3.7)  The  numbers 
in  circles  relate  to  the  nodal-point  numbers.  Those  within 
the  mesh  relate  to  element  numbers.  The  model  is  constructed 
using  ANSYS  STIF42. 

3.5.  Finite  Element  Results 
In  all  the  stress  distribution  results,  the  stresses  are 
normalized  with  respect  to  the  average  applied  axial  stress, 
pav,  acting  on  the  central  adherend.     Accordingly  a,  =  a~/p  , 

£y   "    ^y/Pav/    £xy   =    <*xy/P.v/    While   X   =  X/C. 

where    ax:  axial  normal  stress 

ay:  lateral  normal  stress 
o^:  shear  stress 
pav:  applied  average  stress 
c:  2.56  inches 

For  nodal  stress,  the  stresses  at  each  node  are  computed  by 
averaging  the  nodal  stresses  of  the  node  being  processed  for 
all  elements  that  are  connected  to  that  node,  m  other  words, 
the  component  nodal  stresses  that  are  automatically  stored  per 
element  are  averaged  at  a  node  whenever  two  or  more  elements 
connect  to  the  same  node.  For  example  average  nodal  stress 
ax  at  the  node  N  can  be  calculated  as  below. 


Figure  3 . 6  The  layout  of  finite  element  mesh  in  simple  tension 
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Figure  3.7  The  layout  of  finite  element  mesh  in 
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elementl 

element 2 

-N  1 

element 3 

From  element  1      ax    at  node  N 
From  element  2      ax    at  node  N 
From  element  3      cr_    at  node  N 


average  nodal  stress  =  (alx  +  a2x  +  a3x)/3 

Four  cases  are  considered  in  the  following  four  subsections. 
Case  1  and  2  are  for  two  different  choices  of  adherends  with 
simple  tension  loading.  Case  3  and  4  are  for  the  same  two 
choices  of  adherends  as  case  1  and  2,  but  with  pure  bending 
loading.  Each  of  the  four  cases  is  treated  first  for  double- 
doubler  reinforcement  of  a  continuous  member,  and  then  for  a 
double-doubler  lap  joint. 

3.5.1  Case  1.  Aluminum  Adherends  and  Eooxy  Adhesivg  in 
Simple  Tension 

The  model   for  the   first  case  is  composed  of  aluminum 
adherends    and    an    epoxy    adhesive.         For  double-doubler 
reinforcement  of  a  continuous  member  (Figure  3.1a),  the  axial 
distribution  of  the  lateral  normal  stress  (ay) ,  and  the  shear 
stress   (a^)   along  the  boundary  between  central  adherend  and 
adhesive  are  plotted  from  the  averaged  nodal  stress  values. 
The  reason  why  the  averaged  nodal  stress  values  are  used  is 
that    the    lateral    normal    (ay)    and    shear    stress    {0xy)  are 
continous  across  the  interface.     These  stress  distributions 
are  shown  in  Figure  3.8.      From  these  graphs  and  the  known 
symmetry  about  the  mid-length,   we  can  see  that  the  stress 
distributions    {oy,axy)    along   the   axial    direction   are  very 
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Figure  3.8  Case  1.  The  stress  distributions  of  double-doubler 
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boundary  between  the  central  adherend  an/adhesive 
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uniform  and  at  a  very  low  level  for  -0.75<  x/c  <0.75.  But 
near  the  free  edges  (for  0.75  <  x/c  <  1  and  -1  <  x/c  <-0.75) 
we  can  see  sharp  variation  of  stresses  and  see  the  peak 
stresses  at  the  free  end  of  adhesive.  Different  peak  stress 
values  are  obtained  depending  upon  the  size  of  the  mesh 
element.  This  means  that  calculated  stress  fields  near  this 
area  are  not  quite  reliable. 

In  linear  elastic  analysis,  stress  fields  near  the  two 
bonding  edges  show  singular  behavior.  In  order  to  describe 
the  stress  fields  near  these  corners,  singular  stress 
parameters  are  used  in  Chapter  5. 

The    cohesive    stress    distributions    along    the  axial 
direction   at   different    levels    of   y   of   the    adhesive  were 
obtained.     These  different  levels  of  y  in  the  adhesive  strip 
is  shown  in  Figure  3.9.     According  to  the  results  which  are 
shown  in  Appendix  A,   the  axial  distribution  of  shear  stress 
(0xy)  at  the  level  of  y=o. 06375  shows  a  similar  trend  to  that 
in  the  stress  distribution  along  the  boundary  between  central 
adherend  and  adhesive.    The  axial  distribution  of  axial  normal 
stress    (cx)    at  the  same  level  of  y=o. 06375  shows  a  similar 
trend  as  that  in  the  stress  distribution  along  the  boundary 
between  central  adherend  and  adhesive,  except  that  the  axial 
normal  stress  distribution  drops  toward  zero  near  the  free  end 
of  the  adhesive   layer.      The  axial   distribution  of  lateral 
normal  stress  K)  near  the  free  end  shows  a  similar  trend  to 
that  in  the  stress  distribution  along  the  boundary  between 
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central  adherend  and  adhesive,  except  that  the  peak  stress 
value  becomes  lower  than  previous  peak  stress  on  the  boundary. 

The  axial  distribution  of  axial  normal  stress   (ax)  and 
lateral  normal  stress  (ay)  at  the  higher  level  (y=0.0675)  than 
the  previous  case   (y=0. 06375)   shows  a  similar  trend,  except 
that    the   peak   stress    values    are    lower   than    ones    of  the 
previous  case.     The  axial  distribution  of  shear  stress  (CTxy) 
near  the   free  end   zone,    which   is   taken   from  the   level  of 
y=0.0675,    rises  steeply,    far  beyond  its  average  level,  and 
drops  steeply  near  the  free  end  to  near  zero.     These  stress 
distribution  trends  are  similar  to  the  case  which  is  taken 
from  the  upper  level  of  the  adhesive  at  y=0.08.     The  stress 
distribution  at  the  level  of  y=0.0925  shows  same  trend  as  the 
previous   case   except   that   the   lateral    normal    stress  (CTy) 
changes  its  sign  and  drops  to  a  negative  value  which  means 
compressive  stress,  immediately  after  the  peak  positive  value 
in  the  end  zone. 

The  axial  distributions  of  stresses   (ay,   0xy)   along  the 
boundary  between  adhesive  and  outer  adherend  are  plotted  in 
the  same  way  as  for  the  boundary  between  the  central  adherend 
and  adhesive.     The  results  are  shown  in  Figure  3.10.  They 
show  that  the  shear  stress   (0xy)   drops  to  near  zero  near  the 
free  end.      The  lateral  normal  stress    (ay)    shows  that  after 
the  peak  positive  stress  value  near  the  end  zone,  it  suddenly 
drops  to  the  negative  value  which  means  compressive  stress. 
We  find  that  along  the  boundary  between  the  central  or 
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reinforcement  of  a  continous  member  along  the 
boundary  between  the  outer  adherend  and  adhesive 
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outer  adherend  and  adhesive  the  stress  is  not  well  defined  at 
the  ends  of  adhesive.  Therefore  one  would  not  expect  the  free 
end  condition  (ax=0,  ^=0)  to  be  satisfied  at  the  very  near 
end  zone.  However,  the  end  condition  (o^o^O)  at  the  various 
levels  inside  the  adhesive,  which  means  no  external  traction 
force  in  the  x  or  y  direction,  is  almost  satisfied  even  by 
using  the  displacement  finite  element  method  . 

To   check   the   stress   variation   through   the  thickness 
direction  of  the  adhesive,  the  interlaminar  adhesive  layer  is 
divided  into  four  longitudinal  strips  as  shown  in  Figure  3.9. 
The  results  from  the  data  plotted  for  various  levels  of  y  in 
Appendix  A  show  that  through  the  adhesive  thickness  the  shear, 
axial  and  lateral  normal  stresses  are  pretty  uniform  through 
the  adhesive  thickness    within  the  zone  -0.975  <  x/c  <  0.975, 
but  near  the  free  end  zones  (0.98  <  x/c  <  l  and  -1  <  x/c  <  - 
0.98)     the    stress    changes    rapidly    through    the  adhesive 
thickness.    If  we  assume  that  the  stresses  are  uniform  through 
the  adhesive  thickness,   we  can  not  predict  the  compressive 
lateral  normal  stress  of  the  adhesive  layer. 

For  a  double-doubler  lap  joint  (Figure  3.1b),  the  stress 
distributions  on  the  boundaries  between  central  adherend  and 
adhesive  and  between  outer  adherend  and  adhesive  are  shown  in 
Figure  3.11  and  Figure  3.12,  respectively.  The  axial 
distributions  of  stress  at  different  levels  of  y  are  shown  in 
Appendix  B.  From  the  results  we  find  that  the  axial 
distribution  of  the  lateral  normal  stress  {oy)  ,  axial  normal 
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stress  (ax)  and  the  shear  stress  (CTxy)  are  very  uniform  and  at 
a  very  low  level  between  0.3  <  x/c  <0.8  and  -0.3  <  x/c  <-0.8. 
But  near  the  corner  edge  of  adhesive  line  and  bonded  corner 
at  the  butted  region  between  two  central  adherends,  there  is 
a  very  high  stress  concentration  and  very  complicated  stress 
distribution.    The  peak  stress  values  obtained  from  the  finite 
element  method  are  not  guite  reliable,  since  the  peak  values 
vary  depending  upon  the  mesh  size.    The  stress  distributions, 
except   near  the   bonded   corner   in   the   butted   region,  show 
similiar  trends  to  the  ones  of  double-doubler  reinforcement 
of  a  continous  member.     Near  the  bonded  corner  of  the  butted 
region,    the  shear  stress    (0xy)    is  very  large.      The  lateral 
normal  stress  (ay)  has  its  highest  peak  value  at  the  edge  of 
the  interface  between  central  adherend  and  adhesive.  This 
lateral  normal  stress  (oy)  becomes  compressive  at  the  level  of 
y=0.0925  which  is  very  near  the  interface  between  adhesive 
and  outer  adherend,  and  at  that  interface  it  shows  very  high 
compressive  lateral  normal  stress  (a  )  . 

*    y  / 

3.5.2  Case  2.   Aluminum  Central.  Adhere,   Epoxv  Arih^-i^ 
and  Composite  ggg  Adherent  g  g&gg  < 

The  model  for  the  second  case  is  composed  of  aluminum 
central  adherend,  epoxy  adhesive  and  graph ite/epoxy  composite 
outer  adherend.  A  unidirectional  graphite/epoxy  laminate  is 
used  for  outer  adherend. 
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For  double-doubler  reinforcement  of  a  continous  member, 
the  axial  distribution  of  the  cohesive  lateral  normal 
stress (CTy),  shear  stress (a^)  and  the  axial  normal  stress (ox) 
are  plotted  from  the  averaged  nodal  stress  values  and  are 
shown  in  Appendix  C  for  various  levels  of  y.  The  axial 
distribution  of  the  lateral  normal  stress  and  shear  stress 
along  the  boundary  between  the  central  adherend  and  adhesive 
are  plotted  from  the  averaged  nodal  stress  values  stress  and 
are  shown  in  Figure  3.13.  The  stress  distributions  on  the 
boundary  between  the  adhesive  and  outer  adherend  are  plotted 
in  the  same  way  and  are  shown  in    Figure  3.14. 

For  a  double-doubler  joint,  the  axial  distributions  of 
the  cohesive  stress  at  different  levels  of  y  are  plotted  and 
shown  in  Appendix  D.  The  stress  distributions  on  the  boundary 
between  the  central  adherend  and  adhesive,  and  the  outer 
adherend  and  adhesive,  are  shown  in  Figure  3.15  and  Figure 
3.16,  respectively.  The  stress  distribution  trends  are 
similar  to  Case  1,  but  higher  peak  stress  values  are  obtained 
from  Case  2. 

3'5'3  BendinminUI"  AdherPnds  and  RPQXV  ArihP<n„o  ^ 

The  model  for  Case  3  is  composed  of  aluminum  adherends 
and  epoxy  adhesive  and  is  under  pure  bending  in  the  positive 
sense    (concave   upward).      a   linear  varying   force  boundary 
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Figure  3.16  Case  2.  The  stress  distributions  of  a  double 
doubler  joint  along  the  boundary  between 
the  outer  adherend  and  adhesive 
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condition  equivalent  to  a  couple  is  used  at  the  outer  end  of 
the  central  adherend. 

For  double-doubler  reinforcement  of  a  continous  member, 
the  axial  distributions  of  the  cohesive  stress  are  plotted 
from  the  averaged  nodal  stress  values  and  shown  in  Appendix 
E.    We  have  four  overlap  interfaces  in  the  pure  bending  case. 
The   axial   stress    (ay,    0xy)    distributions   on  the  boundaries 
between  the  central  adherend  and  upper  adhesive  and  between 
the    upper    outer    adherend    and    adhesive    are    plotted  from 
averaged  nodal   stress  values  and  shown   in  Figure   3.17  and 
Figure  3.18,    respectively.      The  axial  stress  distributions 
(ay,   0xy)   on  the  boundaries  between  the  central  adherend  and 
lower    adhesive    and    between    the    lower    outer    adherend  and 
adhesive     are     shown     in     Figure     3.19     and     Figure  3.20, 
respectively.    As  expected,  the  shear  stress  is  symmetric  with 
respect  to  the  mid  plane,  while  the  lateral  normal  stress  is 
approximately  antisymmetric. 

For  a  double-doubler  lap  joint,  the  axial  distributions 
of  the  cohesive  stress  for  various  levels  of  y  in  the  upper 
and  lower  adhesive  are  shown  in  Appendix  F.  The  axial  stress 
distributions  (oy,  »,)  on  the  boundary  between  the  central 
adherend  and  upper  adhesive  and  between  the  upper  outer 
adherend  and  adhesive  are  shown  in  Figure  3.21  and  Figure 
3.22,  respectively.  The  axial  stress  distributions  (0y,  Qxj) 
on    the    boundary    between    the    central    adherend    and  lower 
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adhesive  and  between  the  lower  outer  adherend  and  adhesive  are 
shown  in  Figure  3.23  and  Figure  3.24,  respectively. 
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Figure  3.17  Case  3.  The  stress  distributions  of  double-doubler 
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Figure  3 . 18  Case  3 .  The  stress  distributions  of  double-doubler 
reinforcement  of  a  continous  member  along  the 
boundary  between  the  upper  outer  adherend  and 
adhesive 
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Figure  3.19  Case  3.  The  stress  distributions  of  double-doubler 
reinforcement  of  a  oontinous  member  along  the 
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Figure  3.20  Case  3 .  The  stress  distributions  of  double-doubler 
of  a  continous  member  along  the  boundary 
between  the  lower  outer  adherend  and  adhesive 
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Figure  3.21  Case  3.  The  stress  distributions  of  a  double 
doubler  joint  along  the  boundary  between 
the  central  adherend  and  upper  adhesive 
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Figure  3.24  Case  3.  The  stress  distributions  of  a  double- 
doubler  nomt  along  the  boundary  between 
the  lower  outer  adherend  and  adhesive 
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3.5.4  Case  4.  Aluminum  Central  Adherend .   Epoxv  Adhesive, 
and  Composite  Outer  Adherend  in  Pure  Rending 

The  model  for  Case  4  consists  of  aluminum  central 
adherend,  epoxy  adhesive,  and  graphite/ epoxy  outer  adherend 
under  pure  bending  in  the  positive  sense. 

For  double-doubler  reinforcement  of  a  continous  member, 
the  axial   distributions   of  the   cohesive   stress   at  various 
levels  of  y  in  the  upper  adhesive  are  shown  in  Appendix  G. 
The    axial    stress    distributions    (ay,    av)    on    the  boundary 
between  the  central  adherend  and  upper  adhesive  and  between 
the  upper  outer  adherend  and  adhesive  are  shown  in  Figure  3.25 
and  Figure  3.26,  respectively.    The  axial  stress  distributions 
K,    av)    on  the  boundary  between  the  central   adherend  and 
lower    adhesive    and    between    the    lower    outer    adherend  and 
adhesive     are     shown     in     Figure     3.27     and     Figure  3.28, 
respectively. 

For  a  double-doubler  lap  joint,  the  axial  distributions 
of  the  cohesive  stresses  for  various  levels  of  y  in  the  upper 
adhesive  are  shown  in  Appendix  H.     The  stress  distributions 
K,    av)    on  the  boundary  between  the   central   adherend  and 
upper   adhesive    and   between   the   upper   outer    adherend  and 
adhesive     are     shown     in     Figure     3.29     and     Figure  3.30, 
respectively.       The    stress    distributions    K,     0xy)    on  tne 
boundary  between  the  central  adherend  and  lower  adhesive  and 
between  the  lower  outer  adherend  and  adhesive  are  shown  in 
Figure    3.31    and    Figure    3.32,    respectively.        The  stress 
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distribution  trends  are  similar  to  Case  3,  but  higher  peak 
stress  values  are  obtained  from  Case  4. 
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Figure  3.25  Case  4.  The  stress  distributions  of  double-doubler 
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Figure  3.29  Case  4.  The  stress  distributions  of  a  double 
doubler  joint  along  the  boundary  between 
the  central  adherend  and  upper  adhesive 
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Figure  3.31  Case  4.  The  stress  distributions  of  a  double- 
doubler  joint  along  the  boundary  between 
the  central  adherend  and  lower  adhesive 
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Figure  3.32  Case  4.  The  stress  distributions  of  a  double- 
doubler  joint  along  the  boundary  between 
the  lower  outer  adherend  and  adhesive 


CHAPTER  4 

PARAMETRIC  STUDY  OF  THE  DOUBLE-DOUBLER  JOINT 

4 . 1  Introduction 
The  objective  of  this  chapter  is  to  study  the  parameters 
which  influence  the  stress  distributions,  so  that  designers 
can  design  the  adhesive  bonded  joints  more  efficiently.  The 
ANSYS  finite  element  program  is  used  as  the  computational  tool 
to  determine  the  influence  of  various  joint  parameters  on 
stress  distributions  in  the  double-doubler  joint  under  tesile 
loading. 

The  stress  distributions  presented  in  this  parametric 
study  are  normalized  with  respect  to  the  applied  average 
stress  pav  acting  on  the  central  adherend.  The  influence  of 
the  parameters  is  studied  by  varying  one  parameter  at  a  time 
from  a  double-doubler  joint  which  is  described  in  Figure  3.1. 

The  geometric  and  material  parameters  which  are  examined 

are 

a)  bonding  half  length  c 

b)  the  ratio  of  the  adherend  modulus  to  the 
modulus  of  the  adhesive  E/Ea 

For    the    material    parameter    study,    the    central    and  outer 
adherends  are  made  of  the  same  isotropic  material. 
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The  thickness  of  the  central  adherend,  2tc,  is  constant 
through     the      study.  All      geometric      parameters  are 

nondimensionalized  by  dividing  by  te.  Note  that  x/c  = 
(x/tc)/(c/tc)  .  Material  parameters  are  nondimensionalized  by 
dividing  the  adherend  modulus  E  by  the  modulus  of  the  adhesive 
E.. 

The  influence  of  E/Ea  on  the  stresses  for  simple  tension 
is  dicussed  in  Section  4.2.  The  influence  of  the  length  of 
overlap  on  the  stresses  for  simple  tension  is  examined  in 
Section  4.3. 

4.2  The  Material  ParaTne>1-gT- 
The    material    parameters    consist    of    the  mechanical 
properties   of   the   adherend   and   adhesive.       The  mechanical 
properties   of  the   any   isotropic   elastic  material   are  the 
elastic    modulus    and    the    Poisson's    ratio.       Stresses  are 
affected  by  variation  in  Poisson's  ratio  if  the  changes  in 
Poission's   ratio   are  very   large.      since   the   range   of  the 
variation  of  Poisson's  ratio  for  common  adhesive  systems  is 
very  limited,  only  a  single  material  parameter  is  considered 
in  this  analysis.    This  single  material  parameter  is  the  ratio 
of  the  adherend  modulus  to  the  adhesive  modulus.    To  study  the 
influence  of  the  material  parameters  four  joints  are  analyzed. 
The  adherends  in  the  four  joints  have  four  different  moduli 
of    elasticity    but    the    same    adhesive.        The    modulus  of 
elasticity    is    the    same    for    central    adherend    and  outer 
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adherend.  The  Poisson  ratio  of  the  adherend  is  constant  at 
0.33. 

The  variation  in  the  shear  and  lateral  normal  stress 
distributions  in  the  adhesive  are  examined  for  E/Ea  =  l, 
E/Ea  =  30,  E/Ea  =  loo,  and  E/Ea  =  200.  in  this  study  the 
modulus  of  the  adhesive  Ea  is  0.356xl06  psi  and  the  constant 
Poisson1 s  ratio  of  the  adhesive  is  0.39.  The  stress 
distribution  and  singularity  parameters  in  this  double-doubler 
joint  are  studied  for  four  ratios  of  E/Ea. 

4.2.1  Influence  of  E/Ea 

The  influence  of  the  ratio  of  the  modulus  of  adherends, 
E,  to  the  modulus  of  the  adhesive,  Ea,   on  the  stresses  in  a 
double-doubler  joint  is  examined  here.    The  influence  of  E/Ea 
on  the  axial  distributions  of  shear  stress  at  the  interface 
with  the  center  and  outer  adherend  is  shown  in  Figure  4.2  and 
4.4,   respectively.     Each  figure  consists  of  two  parts.  One 
frame   is   the   neighborhood   of  the   free   end   surface   of  the 
adhesive;  the  other  frame  is  the  neighborhood  of  the  butted 
adhesive  corner.    The  magnitude  of  the  shear  stress  decreases 
rapidly  from  the  neighborhood  of  the  end  of  adhesive  corner 
and  is  almost  uniform,  and  equal  to  the  nominal  shear  stress, 
over  the  middle  three-fifths  of  the  overlap  (0.2  <  x/c  <  0.8 
and  -0.2  <  x/c  <-o.8).     The  shear  stress  then  drops  to  zero 
at  the  free  end  surface  of  adhesive  and  outer  adherend.  As 
E/Ea  increases,  the  shear  stress  concentration  factor  (rmax/rav) 
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decreases  and  the  order  of  stress  singularity  ,A,  decreases. 
(See  chapter  5.3  for  singularity  order  A)     The  magnitude  of 
the  highest  shear  stress  is  greatest  when  E  =  Ea. 
The  location  of  maximum  shear  stress  is  not  affected  when  E/Ea 
is  varied.     The  variation  of  stress  concentration  factor  and 
singularity  order  a  with  respect  to  E/Ea  is  given  in  Table 
4.1.      The  ideal  double-doubler  joint  would  have  a  uniform 
shear   stress   distribution   and  would   also  have   no  lateral 
normal  stress.    The  distributions  of  lateral  normal  stress  on 
the  boundary  between  the  central  adherend  and  adhesive  and 
between  the  adhesive  and  outer  adherend  are  shown  in  Figure 
4.1  and  Figure  4.3,  respectively.    Each  figure  consists  of  two 
parts,    one  frame  is  the  neighborhood  of  the  free  end  surface 
of  the  adhesive;  the  other  frame  is  the  neighborhood  of  the 
butted  adhesive  corner.    The  lateral  normal  stress  is  tensile 
on  the  end  of  boundary  between  central  adherend  and  adhesive 
and  is  compressive  on  the  end  of  boundary  between  adhesive  and 
outer    adherend.       The    lateral    normal    stresses    are  almost 
uniform  and  almost  zero  in  magnitude  over  the  middle  three- 
fifths  of  the  overlap  (-0.2  <  x/c  <-o.9  and  0.2  <  x/c  <  0.9) 
but   rise   rapidly   to   tension   and   compression  peaks   at  the 
respective  ends  of  the  overlap.      The  axial  distribution  of 
lateral  normal  stress  near  the  free  end  of  the  adhesive  layer 
also  have  tension  and  compression  peaks  depending  upon  the 
level  of  y.     it  is  quite  obvious  that  compression  and  tension 
lateral  normal  stress  must  simultaneously  exist  in  order  to 
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satisfy  the  vertical  equilibrium  condition  of  the  adhesive 
layer.     If  we  assume  that  the  stresses  through  the  adhesive 
thickness  are  uniform,    we  can  not  predict  any  compressive 
lateral   normal   stress   in  the  adhesive  layer.      The  tensile 
lateral  normal  stress  concentration  factor  at  the  interface 
between  central  adherend  and  adhesive  is  always  greater  than 
that  at  the  interface  between  adhesive  and  outer  adherend. 
This  concentration  factor  decreases  as  E/Ea  increases.  The 
greatest  concentration  factor  is  obtained  at  both  interfaces 
for  E  =  Ea.     It  is  seen  from  Figures  4.1  to  4.4  that  the 
magnitude   of   the   maximum   lateral    normal    stress    is  always 
greater  than  the  magnitude  of  the  maximum  shear  stress  in  a 
given  joint.     Therefore  the  magnitude  of  the  lateral  normal 
stress  will  be  the  governing  criterion  for  joint  failure  in 
this  double-doubler  joint. 

Table  4.1  Shear  stress  concentration  factor  and  the  order  of 
stress  singularity  a  for  case  l  in  Chapter  5.3 


shear  stress 
concentration  factor 

order  of 
stress  singularity  a 

E/Ea  =  1 

23.2 

0.465 

E/Ea  =30 

3.71 

0.355 

E/Ea  =  100 

1.95 

0.348 

E/Ea  =200 

1.45 

0.346 

respect  to  average  shear  stress 


0.  4 

CD 

W 

1  1  1 

LLf 

cr 

h- 

0.  3 

w 

 i 

< 

cr 

□ 

0  2 

Z 

< 

a: 

LU 
1— 

0.1 

LA' 

Q 

LU 

N 

D  - 

AL 

cr 

o 

-0.1  - 

0 . 

0. 02  -] 

in 

in 

LU 

EX 
1- 

-0.O2  - 

10 

_l 

-0.04  - 

< 

-0 . 06  — 

cr 

O 

-  0 . 08  — 

AL 

-0.  1  - 

cr 

LU 

-0 .  12  — 

< 
_J 

-0. 14  - 

Q 

LU 

N 

-0.16  - 

IAL 

-0.18  - 

;> 

cr 

ON 

-0.2  -- 

1 

1 
1 

1 
) 
( 

 E/Ea=1 

 E/  Ea=30 

/ 

( 

 E/ Ea=lO0 

r 

1 

r  I 

 E/Ea=2Q0 

i  1 
i  1 
/  J 
r  1 

'  1 
<  J1' 
I  /is 

LJ  /i 

r 

'       1        1        1        i  1  r  1  1  1  1  r  

/ 

/ 

s 

0  .  87 


NORMALIZED  AXIAL  COORDINATE  x/c 


0  99 


NORMALIZED  AXIAL  COORDINATE  x/c 
Figure  4.1  Axial  lateral  normal  stress  distribution  on  the 
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Figure  4.2  Axial  shear  stress  distribution  on  the  boundary 
between  the  central  adherend  and  adhesive  for 
four  values  of  E/Ea  or 
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3  Axial  lateral  normal  stress  distribution  on  the 
boundary  between  the  outer  adherend  and  adhesive 
four  values  of  E/E 
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Figure  4.4  Axial  shear  stress  distribution  on  the  boundary 
between  the  outer  adherend  and  adhesive  for  four 
different  values  of  E/E 
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4.3  The  Geometric  Parameters 
The  influence  of  the  geometric  parameters  is  investigated 
by  varying  one  dimension  while  all  others  are  held  the  same 
as  that  of  original  joint.  The  dimensions  of  the  joint  are 
given  in  the  notation  shown  in  Figure  3.2. 

-  half  length  of  overlap:  c  =  2.56 

-  thickness  of  central  adherend:  2tc  =  0.12" 

-  thickness  of  outer  adherend:  tG  =  0.04" 

-  thickness  of  adhesive:  ta  =  0.04" 

The  thickness  of  the  central  adherend  is  held  constant 
throughout  the  analysis.  The  material  parameters  of  the 
adherends  and  adhesive  are  the  same  for  all  the  joints 
analyzed  for  the  study  of  the  influence  of  the  geometric 
parameters : 

-  adherend  modulus:  E  =  10.6E6  psi 

-  adherend  Poisson's  ratio:  v  =  0.33 

-  adhesive  modulus:  E  =  3.6E6  psi 

-  adhesive  Poisson's  ratio  v  =  0.39 

4.3.1  Influence  of  Overlap  T,Pngth 

The  length  of  overlap  is  usually  considered  to  be  the 
most  important  parameter  to  the  designer.  As  the  applied  load 
is  increased,  the  overlap  length  must  be  increased  in  order 
to  reduce  the  magnitudes  of  the  average  stresses.  The  stress 
distributions  in  three  double-doubler  joints  are  analyzed 
here.     All  geometries  of  the  three  joints  are  the  same  as  in 
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original  joint  (c/tc  =  43)  except  the  overlap  lengths  which 
range  between  c/tc  =  30  and  c/te  =  50. 

The  shear  stress  distributions  on  the  boundary  between 
the    central    adherend    and    adhesive    and    between    the  outer 
adherend  and  adhesive  for  the  three  double-doubler  joints 
(c/tc  =  30,   c/tc  =  43,  c/tc  =  50)   show  that  the  shear  stress 
concentration  factors  for  the  three  cases  are  almost  the  same. 
There   is  no   significant  difference   in  stress  distribution 
shape  and  normalized  by  pav  peak  stress  values  on  the  range 
between    c/tc=30    and    c/tc=50.       Although    the  concentration 
factors   are   about  the   same,    the  actual   peak   stresses  are 
significantly  higher   for  the   shorter   overlap,    because  the 
average    stresses    are   higher.       The    lateral    normal  stress 
distribution  on  the  boundary  between  the  central  adherend  and 
adhesive,   and  between  the  outer  adherend  and  adhesive  also 
shows  that  there  is  no  significant  change  on  the  range  between 
c/te  =  30  and  c/tc  =  50.      For  overlap  length  c/tc  =  l  (not 
shown)  the  highest  shear  stress  occurs  at  the  butt  corner  of 
the  overlap  between  two  central  adherends.     As  the  length  of 
overlap  increases  significantly,  the  position  of  the  highest 
shear  stress  approaches  the  ends  of  the  overlap.    The  maximum 
lateral  normal  stress  occurs  at  the  ends  of  the  overlap  for 
all  three  joints. 


CHAPTER  5 
TREATMENT  OF  STRESS  SINGULARITY 
AT  BONDING  CORNERS  IN  THE  DOUBLE-DOUBLER  JOINT 

5.1  Introduction 
In  the  linear  elastic  analysis,   stress  fields  near  the 
two  bonding  edges  between  the  adherend  and  the  adhesive  layer 
show  singular  behavior.     However,   the  finite  element  method 
solutions    provide    finite    valued    stresses    at    the  bonding 
corners  and  in  these  regions  the  maximum  stresses,  calculated 
by  the  finite  element  stress  analysis,  depend  on  the  size  of 
element  meshes  and  are  therefore  not  guite  reliable,    m  order 
to  describe  the  stress  fields  around  these  corners,  singular 
stress  parameters  are  used  to  evaluate  the  magnitude  of  the 
stress  fields  near  the  singularity  point  of  the  double-doubler 
adhesive  joint.      The  assumption  is  that  at  some  small  but 
finite  distance  from  the  singular  point  where  the  material  is 
actually   elastic,    the   singularity   solution   gives   a  closer 
approximiation  to  the  stress  variation  with  distance  from  the 
corner  than  the  finite  element  solution  does. 

This  singularity  solution  method  uses  two  parameters,  ClJ 
and  a  =  i-Pl.  TheSe  parameters  express  the  stress 
distribution  near  a  bonding  edge  along  the  interface,  as 
follows. 
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"H  =  Ci3/rx  5.1 
where 

ai6  is  stress, 

r  is  distance  from  singularity  point, 
0^  is  intensity  of  stress  singularity, 
A  is  order  of  stress  singualrity. 

In  this  chapter  the  ways  to  find  the  and  A  for  double- 
doubler  joints  are  presented.  For  most  practical  purposes  the 
width  of  an  adhesive  joint  is  large  compared  to  the  thickness 
of  the  adherend  and  the  adhesive.  Therefore  plane  strain 
analysis  is  used  throughout. 

5.2  Stress  Singularities  at  Interface  Cnrnprs 
in  Bonded  Dissimilar  Materials 

The  typical  study  of  the  stress  singularity  at  the 
interface  corner  in  bonded  dissimilar  elastic  wedges  has  been 
done  extensively  by  Bogy  [13,   14,  15]. 

Let  D2,  Dl  denote  the  open  two-dimensional  regions 
occupied  by  the  cross  section  of  the  two  wedges  of  angles,  a, 
b,  respectively  (a+b  <  2*) ,  and  suppose  that  they  have  one 
straight  segment  B  of  their  boundaries  in  common  (Figure  5.1)  . 
Denote  their  remaining  straight  boundary  segments  by  Bl,  B2 . 
Let  (r,e)  be  plane  polar  coordinates  of  a  point  with 
rectangular  cartesian  coordinates  (x,y) .  The  material 
properties  are  represented  by  their  respective  shear  moduli 
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y 

A 


Figure  5.1  Two  edge-bonded  elastic  wedges 
of  different  materials  under 
normal  and  shear  loading 
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and  Poission's  ratios  (G^i/i)  and  (G2,i/2).  Shear  and  normal 
tractions  are  applied  to  their  free  boundaries:  tx(r)  and 
nx(r)  on  Bl,  and  t2(r)  and  n2(r)  on  B2.  The  traction  is 
independent  of  the  coordinate  perpendicular  to  the  cross 
section.  Therefore  we  consider  this  boundary  value  problem 
in  the  plane  theory  of  linear  elastostatics.  The  problem  now 
is  to  find  the  stress  and  displacement  fields  with  plane  polar 
components  related  to  the  Airy  stress  function  $X/  $2,  which 
are  suitably  defined  and  satisfy  the  bi-harmonic  equation 


V*  *  =  0  in  Dl,   D2  5>2 


The  stresses  at  any  point  (r,e)  can  be  derived  by 


l  a<f>  i  a2</> 
 +     


r     ar  r2  ae2 

a2^ 


2  , 


dr 


1        d2<(>  l  d<f) 


°r9    "   + 


r       arae         r2  ae 

aur  i        i  d*       i  d24 

=  (         +  — ,  — :   -(i  -  m/4)vV) 

dr  2G       r  ar       r2  ae2 

where  m  =  4(1-!/)   for  plane  strain 
=  4/(i-i/)  for  plane  stress 


5.3 
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aue        ug      i  aur       i  i  d2<f>        i  a<t> 

3r  r       r  39       g  r  arae       r2  ae 


which  satisfy  following  boundary  conditions  on  Bl,  B2 
(*i)ee(r,-b)  =  nx(r)  (^2)ee(r,a)  =  n2(r) 

(^i)re(r,-b)  =  tx(r)  (a2)re(r,a)  =  t2(r) 


5.4 


as  well  as  the  conditions  of  bonding  on  B  (continuity  of 
traction  and  displacement  components) 

(*i)ee(r,0)    =    (a2)6e(rf0)      (aJ^f^O)    =  (a2)r6(r,0) 

5.5 

(ujr(r,0)  =  (u2)r(r,0)         (uJ.fr^O)  =  (u2)e(r,0) 

In  these  equations  the  subscripts  1  and  2  indicate  quantities 
in  wedges  I  and  II,  respectively.  m  addition  we  shall 
require  the  stress  components  of  the  two  stress  fields  to 
specify  the  regularity  conditions 


°rr,  aze  -  0(r~    )  as  r  t  »    for  every  s  >  0  55 
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The  Mellin  transform  is  the  most  convenient  tool  for 
solving   the    boundary-value   problem.    We    denote    the  Mellin 
transform  of  a  function  f  defined  and  suitably  regular  on 
(0,  «)  by 

g(f'*S)  =  L  f  (r)   r*  ldr  5'7 
where  s  is  the  complex  transform  parameter. 

Let   t(a,B),    S„(s,e),    aM(afe)f    are(s,e),    ur(s,e),  u.(s,e)( 
n(s),     €(s)   in  this  order  denote  the  Mellin  transforms  with 
respect   to   r   of    *(r,e),    r2arr(r,e),    r^^e),  r2ar9(r,e), 
rur(r,e),  rue(r,e),  r2n(r)  ,  and  r2t(r). 

A  formal  application  of  the  Mellin  transform  to 
Equation  5.2  yields  an  ordiary  differential  equation 

d2  d2 
.2 


+  s    3  C  —T    +  (s  +  2)2  ]    $(s,e)  =  0  5.8 


The  general  solution  for  this  differential  equation  is  given 
by  Bogy  [13,  14,  15]  by 


$(s,e)     =  a(s)sin(se)  +  b(s)cos(se)  +  c(s)sin(s  +  2 


e) 


+  d(s)cos(s  +26)  5  9 
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The  functions  a(s)  ,  b(s)  ,  etc.,  (al,  bl,  etc.,  for  £  and 
a2,  b2,  etc.,  for  £2)  are  to  be  determined  through  the 
transforms  of  Equation  5.3  from  the  transforms  of  the  boundary 
condititons  Equation  5.4  and  Equation  5.5.  After  use  of 
Equation  5.7  these  transformed  equations  appear  as 

a  a2 

*rr(S,e)     =     (       -fQ2  ~    S)  £(S,0) 

£ee(s,6)  =  s(s  +  1)  $(s,e) 

$«e(s,e)  =  (s+1)  -|e$(sre)  5.10 
ur(s,e)  =  1/2G  [s  £(s,e)  +mc(s)sin(se+2e)+md(s)cos(se+2e) ] 

ue(s,e)  =  1/2G  [-  _d_  $(s,e)  +  mc(s)cos(se  +  29)  - 

md(s) sin(se+26) ] 


(  *i)ee(s,-b)   =  njs)  (  Sj^s^b)  =  tl(s) 

5.11 

(  S2)ee(s, a)  =  n2(s)  (  S2)re(S,a)  =  €2(S) 


(  ^i)ee(s,0) 
(  ^,(8,0) 


"  (  ^2)ee(s,0) 
=   (  u2)r(s,0) 


(  *l)r9(S,0) 

(  U2)e(s,0) 


=     (  *2)re(S,0) 

5.12 

=   (  Ui)8(s,0) 
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By  substitution  of  Equation  5.9  and  Equation  5.10  into 
Equation  5.11  and  Equation  5.12,  a  system  of  eight  equations 
for  the  eight  unknown  functions  al(s),  bl(s),   cl(s),  dl(s), 
a2(s),    b2(s),    c2(s),    d2(s)     is    obtained.       Then    Bogy  is 
concerned  with  determining  the  dependence  of  the  order  of  the 
singularity  in  the  stress  fields  as  r  -»  0  on  the  angle,  a,  b 
and  the  material  properties  Glr   ux,   G2,   u2.      It  is  shown  in 
reference     [13,     14,     15]     that    the    stress    fields    have  a 
singularity  of  order  r'^log  r,  or  1  as  r  -»  0      according  as 
the  determinant  of  the  coefficient  in  the  system  of  eight 
equations  has  a  zero  s  in  -2  <  Re(s)<  -l,  has  no  zero  in  this 
open  strip  but  s  =  -2  is  a  zero  of  order  two,  has  no  zero  in 
the   strip   and   s   =   -2    is   only   a   simple    zero.  Therefore 
attention  is  focused  on  the  determinant  of  the  coefficients 
in  the  system  of  equations  with  a  view  toward  locating  its 
zeros  in  -2  <  Re(s)   <  -i.     By  letting  Q  denote  the  value  of 
the  determinant  of  the  coefficients  of  these  eight  equations 
the  following  equation  is  obtained. 

Q(a,b,a,fl;p)  =  A(a,b;p)B2  +  2B(a,b;p)aB  +  C(a,b;p)a2 

+2D(a,b;p)B  +  2E(a,b;p)Q  +  F(a,b;p)  5.13 

where 

A(a,b;p)  =  4M(p,a)M(p,b) 

B(a,b;p)  =  2p2sin2(a)M(p,b)  +  2p2sin2(b)M(p,a) 
C(a,b;p)  =  4p2(p2-i)sin2(a)sin2(b)  +  M[p,  (a-b)  ]       5. 14 
D(a,b;p)  =  2p2[sin2(a)sin2(pb)  -  sin2 (b) sin2 (pa)  ] 


97 

E(a,b;p)  =  -D(a,b;p)  +  M(p,b)  -  M(p,a) 
F(a,b;p)  =  M[p, (a+b) ] 

in  which  the  auxiliary  function  M(p,x)  is  defined  by 

M(p,x)  =  sin2(px)  -  p2sin2(x)  5.15 

In  Equation  5.13,  Equation  5.14  and  Equation  5.15 
p  =  p  +  iri  is  a  complex  varible. 


a  = 


Gim?  "  G?mi  a  _   G,(m,  -  2)   -  G,fm,  -  2) 

+  '       ~  Glm2  +  G2m,   5'16 


where 

in  -  4(1  -  v)  for  plane  strain 

=  4/(1  +  u)   for  generalized  plane  stress 


If  we  let  Q=0  (Equation  5.13)  and  solve  for  the  smallest  value 
of  p,  which  is  called  singularity  parameter,  Pl  ,  then  the 
orders  of  singularity  in  the  stress  field  of  wedges  I  and  II 
are  found  to  be 

°i,  o2  =    0(rp  l)  if  p=Pl  is  real 

Otr^cosf^log  r)  or  r^sin^log  r)  ]  5.17 
lf  Pi  =  £1  +  »?i  is  complex 

The  value  of  the  singularity  parameter  Pl  in  Equation  5.17  is 
dependent  upon  the  two  composite  material  constants  *  and  6 
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from  Equation  5.16.  Stress  near  the  corner  point  can  thus  be 
expressed  for  real  px=p  as 

°ti  =  c"ijfij(e)rp-1  5.18 
where  c"u  is  a  constant  and  f^fe)   is  a  function  of  e. 
Along  the  interface,  e  =  0,  stresses  do  not  depend  on  e,  and 
we  can  write 


"  c  ijrP  1  5.19 
where  c'y  =  C  ^(0)  is  a  constant.    In  component  form  for 
e  =  0    Equation  5.19  can  be  written  as  follows  for  real  Pl=p 

ox  =  arr  =  c»xrp-1 


5.20 


o„  -  a  re  =  C'  rp_1 


5.3  Calculation  of  Order  and  intensifies  of  Stress 
Singularities  at  the  Bonding  Corners  nf  Doubl e-noubler  Jninr 

In  a  double-doubler  joint  model  we  can  specialize  the 
foregoing  discussions  for  the  two  different  cases. 
5.3.1  Case  1.   Half  Plane  Rpnded  to  a  Quarts  Piano 

We  have  four  different  locations  for  this  case 
(Figure  5.2).     However,  the  singularity  parameter  Pl  has  the 
same   value   at   each   of   the    four   locations    for   the  simple 
tension  case  because  of  the  symmetric  loading  and  geometric 
conditions. 

The  expressions  in  Equation  5.14  can  be  written  as  follows: 
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A  =  4sin2(pw)  [sin2(pw/2)  -  p2] 
B  =  2p2sin2(pjr) 

C  =  sin2(p*/2)  5>21 
D  =  -2p2sin2(p*) 

E  =  (2p2  -  l)sin2(p*)  +  sin2(p,r/2)  -  p2 
F  =  sin2(3pw/2)  -  p2 


Substitute  this  Equation  5.21  into  the  Equation  5.13.  Then 
Equation  5.13  becomes 


Q(*,>r/2,a,B;p)  =  4sin2(pw)  [sin2(p7r/2  -  p2]B2 

+  2x2p2sin2(p7r)a6  +  sin2(p^/2)Q2  5.22 
+  2[-2p2sin2(P}r)]B  +  2[(2p2  -  l)sin2(p,r) 


+  sin2(pw/2)   -  p2]  +  sin2(3pJr/2)   -  p2 


The  solution  of  Q=0  in  Equation  5.13  for  the  smallest  p 
yields  the  singularity  parameter  Pl.  This  parameter  Pl  depends 
on  only  two  factors;  the  material  constants  a  and  B.  From  the 
given  material  properties  (Table  5.1)  ,  we  find  the  value  «  and 
B  using  the  Equation  5.16. 

Table. 5.1.  Selected  material  properties  for 


material 

Young's  Modulus 

Poission's  ratio 

Aluminum 

10.6E+6  psi 

0.33 

Epoxy 

3.56E+5  psi 

0.39 
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A  computer  program  was  written  to  solve  the 
characteristic  equation  (Equation  5.22).  The  Muller  method 
[42]  was  used  for  this  computer  program.  This  method  is  used 
to  find  any  prescribed  number  of  zeros,  real  or  complex,  of 
an  arbitrary  function. 

The    values    found    for    a,     B    are    -0.93     and  -0.165, 
respectively.     Using  these  a  and  6  values  we  find  that  p1  = 
0.645.     Then  Equation  5.20  can  be  rewritten  in  the  following 
form  by  dividing  by  the  average  applied  axial  normal  stress 
(Pav  =  P/2tc)  acting  on  the  central  adherend.    In  the  case  that 
r  goes  in  the  negative  x-direction  along  the  interfacial  plane 
between  central  adherend  and  adhesive  (r1  direction  in  Figure 
5.2),  the  Equation  5.20  can  be  written  in  the  follwing  form. 
Here  the  axial  normal  stress  ax  is  not  considered  since  it  is 
not  well  defined  on  the  interface  boundary  between  central 
adherend  and  adhesive. 

,  °y    .  r      Pi-l  r  -0.355 

(      —      )  =  Cy       (    )  =Cy(    ) 

Pav  C  c 

5.23 

,     Qxy  .                       r      Pi-l                      r  -0.355 
(  —  )     "  (    )  =    Cxy  (    ) 

Pav  C  c 

For  the  case  that  r  goes  in  the  positive  y  direction  at  point 
B  (r2  direction  in  Figure  5.2),  the  Equation  5.20  can  be 
written  in  the  following  form.  Here  we  assume  that  joint 
boundary  BB-  is  parallel  to  the  y  axis  at  x=2.56.    It  is  found 
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that  ax,  ff^y  are  almost  zero  near  the  free  end  of  adhesive 
layer  since  there  is  no  traction  force  in  the  x  or  y  direction 
on  this  free  end  of  the  adesive  layer.  Therefore  the 
intensity  of  stress  singularities  Cx  and  are  very  small. 
Only  the  ay  field  is  considered. 

°y  r      p^l  r  -0.355 

(  —  )=Cy(    )  =Cy(    )  5.24 

P.v  ta  ta 

cx/    Cy/    Cxy    are    called    the       "intensities    of  stress 
singularity",  which  are  constant  for  a  specific  joint  under 
a  specific  loading  conditions.     If  we  know  these  factors  for 
a  joint,   the  stresses  near  the  bonding  corner  along  a  bond 
can  easily  be  found  by  using  Equation  5.23  and  Equation  5.24. 
These  explicit  equations  give  infinite  stress  right  at  the 
corner.      The  stress  intensity  parameters  can  be  used  as  a 
means    of    failure    prediction    by    comparing    the  calculated 
parameters  with  known  critical  values  of  these  parameters  for 
the  material  and  comparison  of  stresses  very  near  the  bonding 
corners  where  stress  singularity  occur.     In  order  to  find  the 
desired   intensities   of   stress    singulaity   for   the  double- 
doubler  joint,  the  following  procedure  is  used. 

1)  Using  the  specified  boundary  displacement  (SBD)  method 
described  in  the  following  section  ,  calculate  stress  values 
for  ax,  a^,  ay  from  the  end  of  the  adhesive  interface  corner 
(r/c  =  0)   to  several  nodal  points  in  the  rlt   r2  directions, 
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which  are  very  near  the  end  of  the  adhesive  interface  corner. 
(See  Figure  5.2) 

2)  Substitute  these  values  into  equations  Equation  5.23  and 
Equation  5.24,  and  find  several  values  for  Cx,  Cy,  C^, 
respectively. 

3)  Choose  the  maximum  values  from  the  computed  values  for  Cx/ 

The  maximum  values  are  taken  as  the  desired  intensities  of 
the  stress  singularity.  Substitution  of  these  values  into 
Equation  5.23  and  Equation  5.24  gives  a  set  of  stress 
distributions  near  the  singular  point. 

5.3.2  Case  2.   Three  Quarter  Plane  Bonded  to  a  Quarter  Plang 
(composite  full  plane  a=3w/2,  b=jr/2) 

In  order  to  apply  the  same  general  solution,  Equation 

5.13,  to  the  neighbourhood  of  the  butted  end,   it  is  assumed 

that  the  end  of  the  central  adherend  is  de-bonded  from  the 

adhesive    between    the    two    central    adherends.        This  is 

equivalent  to  assuming  a  straight  crack  that  meets  the  plane 

of  bond  at  a=3,r/2  and  b=*/2.     We  then  have  four  different 

singularity  locations  as  with  Case  1.    (See  Figure  5.3) 

The  singularity  parameter  Pl  has  the  same  value  at  each  of  the 

four  locations   for  the  simple  tension  case  because  of  the 

symmetric  loading  and  geometric  conditions  as  Case  l.  The 

expressions  in  Equation  5.14  can  be  written  as  follows. 
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A  =  4   [sin2(3p*/2)    -  p2]  [sin2(p,r/2)-p2] 
B  =  2p2[sin2(p*/2)  +  sin2(3pw/2)  -2p2] 
C  =  4p2(p2-l)    +  [sin2(P}r)] 

D  =  2p2[sin2(p*/2)  -  sin2(3pw/2)]  5.25 
E  =  (2p2-l)  [sin2(3p*/2)   -  sin2(p*/2)] 
F  =  sin2(2pw) 

Substitute  this  Equation  5.24  into  the  Equation  5.13.  Then 
Equation  5.13  becomes 

Q(3*/2,,r/2,ar,B;p)   =  4   [sin2(3p,r/2)  -p2]  [sin2(pw/2)  -p2]B2 

+  4p2[sin2(pJr/2)+sin2(3p*/2)+2p2]aB  5.26 

+  [4p2(p2+l)  +  sin2(pff)]Q2 

+  4p2[sin2(pff/2)  -  sin2(3p*/2)]B 

+   (4p2+2)  [sin2(3p*/2)   -  sin2(pw/2)]Q 

+  sin2(2p7r) 


The    solution    of    Q    =    0    for    the    smallest    p    yields  the 
singularity  parameter  Pl.     The  values  found  for  a,  B  are 
a  =  -0.93  B  =  -0.165. 

Using  these  a  and  6  values,  we  find  from  Equation  5.26 
that  Pl  =  0.543.  The  same  SBD  procedure  is  used  to  find  the 
desired  intensities  of  stress  singularity  as  was  used  for  Case 
1.  Then  Equation  5.20  can  be  rewritten  in  the  following  form 
in  same  way  as  Case  l.  For  the  case  that  r  goes  along  the 
positive  x  direction  of  the  interface  between  central  adherend 
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and  adhesive  from  B"  (r3  direction  in  Figure  5.3),  Equation 
5.20  can  be  written  in  the  following  form.  Here  only  the 
axial  shear  stress  is  considered  since  it  is  the  dominant 
stress  field. 


°xy  r      px-i  r  -0.457 

  )     =  Cxy(    )  =  Cxy   (    )  5.27 


For  the  case  that  r  goes  along  the  positive  y  direction  from 
B"  (r«  direction  in  Figure  5.3),  Equation  5.20  is  rewritten  in 
the  following  form.  Axial  normal  stress  (ax),  shear  stress 
(ct^)  are  considered  here  since  these  are  the  major  stress 
fields  here. 


,  ax  r      px-i  r  -0.457 

(  —    )  =  Cx  (  —  )  -  C,  (    —  ) 

Pav  t.  t. 


5.28 


ff*y  ,  r  Pi-1  r  -0.457 
—  )      =  C*y    (  —  )  =  (    ) 

Pav  t.  t. 


5.4  Specified  Boundary  pjspl  cement  Mat-hod  rSRm 
5.4.1  Introduct  i  on 

In  many  applications  of  finite  element  structural 
analysis,  it  becomes  necessary  to  determine  concentrated 
stress  in  a  small  portion  of  a  structure.  The  stress  field 
surrounding    a    sharp    fillet    is    characterized    by  rapidly 
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changing  stress  gradients,  which  call  for  a  comparatively  fine 
mesh.  It  can  be  difficult  to  include  sufficient  mesh 
refinement  in  an  entire  model  for  accurate  stress  computation 
at  the  concentration.  One  of  approaches  to  this  problem  is 
the  specified  boundary  displacement  method  (SBD) . 

The  specified  boundary  displacement  method  is  a  technique 
whereby  the  boundaries  of  a  specialized  model  of  a  sub-region 
within    a    structure    assume   the    displacement    field    from  a 
previous  solution  of  the  entire  structure.     The  analyses  of 
the  subregion  and  of  the  entire  structure  are  separate.  What 
is    important   is   that   the    initial    analysis   of   the  entire 
structure  yields   an  accurate  displacement   field  where  the 
boundaries  of  the  subregion  occur  and  that  the  analysis  of  the 
subregion  accurately  recovers  stress  and  strain  at  the  points 
of    interest.       Analyses    of    structures    containing  stress 
concentration  are  excellent  applications   for  the  specified 
boundary  displacement  method.      In  reality,    it  can  be  often 
observed  that  the  more  localized  is  the  effect  of  a  stress 
concentration,  the  greater  is  the  need  for  a  specialized  model 
to  recover  accurate  stress  and  strain.    See,  for  example,  [25, 
26,  33]. 

5.4.2  Procedure 

Analysis  of  a  stress  concentration  by  the  SBD  method 
proceeds  as  described  below  and  shown  in  Figure  5.4.  This 
procedure  follows  that  in  [33]. 


spec  i  f  i  ed 


boundar  ies 


Figure  5.4  Specified  boundary  displacement  method 
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1)  A  finite  element  model  of  a  large  portion  of  the 
structure  is  prepared.  This  model  must  have  sufficient 
refinement  to  obtain  the  overall  displcement  field,  but 
need  not  have  sufficient  refinement  for  accurate 
computation  of  stress  at  the  concentration.  A  stress 
solution  is  performed  for  the  desired  loading  with  this 
coarse  model. 

2)  A  detailed  finite  element  model  of  the  stress 
concentration  is  prepared.  This  submodel  extends  some 
arbitrary  distance  away  from  the  concentration. 

3)  Displacements  along  the  specified  boundaries  of  the 
submodel  are  obtained  from  the  coarse  model  solution.  For 
each  node  along  the  specified  boundary,  the  appropriate 
element  of  the  coarse  model  is  determined.  The  corner 
displacements  of  that  element  are  interpolated  using  its 
shape  function  to  produce  the  displacement  constraints 
to  be  applied  to  the  specified  boundary  nodes. 

4)  The  submodel  stress  solution  is  performed  using  these 
interpolated  displacements  along  with  any  other 
appropriate  boundary  conditions. 

The  SBD  method  permit  independent  definition  of  the 
coarse  and  sub  models.  it  allows  review  of  coarse  model 
results  before  decisions  are  made  about  local  models. 
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5.4,3  Benefits  of  Submodel ling 

Acoording  to  [25,26],  the  benefits  of  submodelling  are 

that 

1)  It  reduces  the  need  for  complicated  transition  regions 
in  solid  finite  element  models. 

2)  It  enables  the  analyst  to  decide  which  areas  of  a 
structure  will  receive  detailed  analysis  after  the 
initial  analysis  is  complete. 

3)  It  enables  the  analyst  to  study  the  effect  of  local 
geometric  changes  around  the  stress  discontinuity  of 
alternate  designs. 

4)  It  enables  the  analyst  to  demonstrate  the  adequacy  of 
mesh  refinement. 

5)  The  refined  mesh  boundary  doesn't  have  to  follow 
coarse  mesh  element  boundaries. 

5.4.4  Coarse  Model  and  Submodel 

Accurate   results   can   be   obtained   through   the  proper 

selection  of  coarse  model  and  submodel.  Two  different  coarse 
models  and  three  different  sub-models  were  prepared  for  the 
double-doubler  joint.  From  these  models  the  combination  of 
coarse  model  and  submodel  that  gave  the  highest  peak  stresses 
was  selected.  Two  cases  were  considered.  The  submodel  for 
Case  1  is  in  the  corner  at  the  free  end  of  the  adhesive  as 
shown  in  Figure  5.5.  For  Case  2  the  submodel  is  in  the 
adhesive  near  the  butted  corner  of  the  central  adherend,  as 
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shown  in  Figure  5.6.     In  both  cases  the  joint  is  loaded  in 
simple  tension.     The  selected  coarse  and  submodel  for  Case  1 
are  shown  in  Figure  5.5.    For  Case  2  the  same  coarse  model  is 
used  as  for  Case  1  and  another  submodel  is  generated  as  shown 
in  Figure  5.6.    The  analysis  was  carried  out  using  the  finite 
element  computer  program  ANSYS.    For  the  model  which  is  under 
simple  tension  loading,  only  one  guarter  of  the  double-doubler 
joint  is  selected  for  this  finite  element  model  because  of 
symmetry.    For  the  model  which  is  under  pure  bending,  half  of 
the  double-doubler  joint  will   be   selected   for  this  finite 
element    model.         Both    coarse    model     and     sub-model  are 
constructed     using  ANSYS  STIF42,   a  four-noded  isoparametric 
plane  strain  element.     Using  the  stress  values  obtained  from 
the  sub-models  and  substituting  these  values  into  Eguation 
5.23,    Eguation   5.24,    Eguation   5.27    and   Eguation   5.28,  we 
obtain  several  different  values  for  Cx,  Cxy/  Cy  for  Case  1  and 
Case  2  for  two  different  loading  cases,   respectively.  Here 
only  the  simple  tension  loading  case  is  presented  to  obtain 
the  values  for  Cx,    C^,    C,  because  the  procedure  to  obtain 
these  values  for  pure  bending  case  is  exactly  the  same  as  for 
simple  tension.     The  obtained    different  values  are  shown  in 
Table  5.2  and  Table  5.3  for  Case  1,   and  in  Table  5.4,  and 
Table  5.5  for  Case  2  for  the  simple  tension  loading. 
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Figure  5.5  The  Coarse  model  and  sub-model 
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Table  5.2.  Different  values  Cy,         for  Case  1  in  rx 
direction  in  simple  tension  loading  case 


*1 

<VPav 

CTxy/Pav 

.002 

0.00079 

0.2066 

0.084 

0.014 

0.0057 

.004 

0.00157 

0.1434 

0.071 

0.013 

0.0065 

.006 

0.00236 

0.1132 

0.068 

0.012 

0.0072 

.008 

0.00314 

0.09 

0.065 

0.01 

0.0077 

.01 

0. 00394 

0.074 

0.06 

0.0095 

0.0077 

S?6"^1  ^S  the  distance  originating  from  the  end  corner  of 
the  adhesive  to  the  negative  x  direction  of  interface  plane 
(See  Figure  5.2)  v 


Table  5.3.  Different  values  Cx,  Cyl  Cjcy  for  Case  2  in  the 
r2  direction  in  simple  tension  loading  case 


r2 

tfy/Pav 

cx 

.002 

0.033 

0.02 

0.3312 

0.0057 

0.094 

.004 

0.067 

0.018 

0.1884 

0. 0066 

0.069 

.006 

0.1 

0.002 

0.138 

0.0008 

0.059 

.008 

0.133 

0.  002 

0.108 

0.0009 

0.051 

.01 

0.167 

0 

0.084 

0.0 

0.043 

nf  fho  •      ti1ST:anc.e  originating  from  the  end  corner 

of  the  adhesive  to  positive  y  direction  (See  Figure  5.2) 
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Table  5.4.  Different  values        for  Case  1  in  the  r3 
direction  in  simple  tension  loading  case 


r3 

r3/c 

''xy/Pav 

.004 

0.00156 

0.273 

0.014 

.008 

0.00312 

0.22 

0.0157 

.012 

0.00469 

0.178 

0.0153 

.016 

0.00625 

0.156 

0.0153 

.02 

0.0078 

0.141 

0.0153 

where  r3  is  the  distance  originating  from  the 
butted  corner  point  to  the  positive  x  direction 
of  interface  plane  (See  Figure  5.3) 


Table  5.5.  Different  values  Cx,         for  Case  2  in  the  r4 
direction  in  simple  tension  loading  case 


r. 

r,/t. 

^x/Pav 

^xy/Pav 

cx 

.0025 

0.0625 

0.3612 

0.1512 

0.102 

0.04 

.005 

0.125 

0.2712 

0.07 

0.1048 

0.027 

.01 

0.25 

0. 1888 

0.046 

0.1 

0.024 

.015 

0.375 

0. 1332 

0.038 

0.  08 

0.  024 

a  *£  tne  ^stance  originating  from  the  butt, 
point  to  the  positive  y  direction  (See  Figure  5.3) 
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From  these  values,  we  choose  the  maximum  values  in  each 
direction.  The  maximum  values  for  Case  1  in  the  r1  direction 
are  Cy  =  0.014,  C^  =  0.0077  and  in  the  r2  direction  Cy  =  0.067. 
For  Case  2  the  maximum  values  in  the  r3  direction  are  C  = 
0.016  and  Cx  =  0.105,  and  in  the  r4  direction  C^  =  0.04. 
Substituting  these  values  into  Equation  5.23,  Equation  5.24, 
Equation  5.27  and  Equation  5.28,  we  obtain  the  final 
equations,  which  give  a  set  of  stress  distributions  near  the 
singular  point  as  follows. 

Case  1  for  the  simple  tension  case  in  the  direction 

-0.355 

(<VPav)    =    0.014  (Ti/C) 

5.29 

(**y/P.v)   =  0.0077  (ri/c)"0-355 

Case  1  for  the  simple  tension  case  in  the  r2  direction 

K/Pav)  =  0.094  (r^tj"0-355  5  3Q 

Case  2  for  the  simple  tension  case  in  the  r3  direction 

("xy/P«v)  =  0.0157(r3/c)-°*57  5  <31 

Case  2  for  the  simple  tension  case  in  the  r,  direction 
(*x/Pav)  =  0.105  (r,/ta)-°457 

5.32 

(**y/Pav)    =0.04  (r4/ta)-°*57 


CHAPTER  6 
RESULTS  AND  DISCUSSION 


6  . 1  Introdntrt-  i  nn 
The    results    of    finite    element    stress    analysis    of  a 
double-doubler  adhesive  bonded  joint  are  presented  in  this 
thesis.     Two  different  loading  cases  (a  simple  tension  and  a 
pure    bending    applied    at    the    ends    of    the    adherends)  are 
considered    to   evaluate   the   traction    distributions    on  the 
interfaces  and  the  cohesive  stresses  in  the  interior  of  the 
adhesive  layer.     Results  of  a  parametric  study  of  the  joint 
properties  is  presented  to  show  the  influences  of  the  various 
joint  parameters  on  the  stress  distributions  in  the  joint.  A 
refined   treatment   for  the   stress    field   near  the  singular 
points  is  presented.  A  set  of  "stress  singularity  parameters- 
are  devised  to  represent  the  stresses  in  the  regions  close  to 
and  including  the  singular  points. 

The  summary  of  the  results  of  this  analysis  is  presented 
in  Section  6.2.  Suggestions  for  more  efficient  joint  design 
are  discussed. 

Recommendations  for  future  work  are  in  the  concluding 
section,  Section  6.3. 
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6.2  Summary  and  Conclusions 
The  double-doubler  joint  with  identical  adherends  and 
non-identical  adherends,   shown  in  Figure  3.1,  was  analyzed. 
The  computational  tool  for  the  analysis  is  the  ANSYS  finite 
element  program.    All  three  materials  involved  are  assumed  in 
linear  elastic  under  the  applied  load.     The  stress  analysis 
is  treated  as  a  plane  strain  problem.  A  linear  elastic  stress 
analysis  of  the  double-doubler  joint  is  given  in  Chapter  3. 
A  parametric  study  of  the  double-doubler  joint  is  given  in 
Chapter    4,    where    the    influences    of    two    important  joint 
parameters  are  analyzed.     In  Chapter  5,   "stress  singularity 
parameters"  are  devised  to  describe  stress  fields  near  the 
singularity  zone. 

The  important  stresses  in  a  double-doubler  joint  are  the 
shear  and  lateral  normal  stresses  in  the  adhesive  and  axial 
normal  stresses  in  the  adherends.  The  maximum  magnitudes  of 
these  stresses  and  their  gradients  occur  at  the  end  zone  of 
adhesive  overlap  and  at  the  bonded  corner  zone  at  butted 
region  between  two  central  adherends  in  the  simple  tension 
case.  In  the  pure  bending  case  the  maximum  magnitude  of  these 
stresses  and  their  gradients  occur  at  the  end  zone  of  adhesive 
over  lap.  The  critical  areas  of  the  adhesive  are,  therefore, 
at  the  ends  of  the  overlap  and  at  the  bonded  corner  of  butted 
region. 

These  critical  areas  extend  over  approximately  10%  of  the 
overlap   length  at  the  ends   of  the   overlap  and   10%   of  the 
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overlap  length  at  the  bonded  corner  of  the  butted  region.  The 
axial  distributions  of  stress  (axf  CTy,  a^)  in  the  adhesive  are 
almost  uniform  at  a  low  level  over  the  middle  75%  of  overlap. 
The  axial  normal  stress  distribution  in  the  adhesive  layer 
varies  linearly  through  the  adhesive  layer  thickness  near  the 
end.  Along  the  boundary  between  the  central  adherend  or  outer 
adherend  and  adhesive,  axial  normal  stress  (aj  is  not  well 
defined.  Therefore  one  would  not  expect  the  free  end 
condition  (ax  =  0xy  =  0)  to  be  satisfied  at  the  very  near  end 
zone.  However,  the  end  condition  (ax  =  =  0)  at  the  various 
levels  inside  the  adhesive,  which  means  no  external  traction 
force  in  the  x  or  y  direction  on  the  end,  is  almost  satisfied 
by  using  the  displacement  finite  element  method. 

The  lateral  normal  stress  has  its  highest  peak  value  at 
the  edge  of  the  interface  between  the  central  adherend  and 
adhesive.  The  lateral  normal  stress  (ay)  becomes  compressive 
in  part  of  the  upper  level  of  adhesive  near  the  interface 
between  the  adhesive  and  outer  adherend,  and  at  that  interface 
it  shows  an  interval  of  very  high  compressive  stress,  if  we 
assume  that  stresses  through  the  adhesive  thickness  are 
uniform,  we  can  not  predict  the  compressive  lateral  normal 
stress  regions  of  the  adhesive  layer. 

From  the  parametric  study  of  the  joint  with  adherends  of 
the  same  isotropic  material,  it  is  concluded  that  the  single 
most  important  parameter  influencing  the  stresses  is  the  ratio 
of  the  adherend  modulus  to  the  modulus  of  the  adhesive,  E/E  . 
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The  stresses  in  the  adherends  and  adhesive  are  influenced  by 
change  in  E/Ea.  The  stress  peaks  are  greatest  when  E/Ea  =  l 
and  decrease  with  increase  in  E/Ea.  The  stresses  in  the 
adhesive  and  adherends  do  not  change  appreciably  for  values 
of  E/E,  >  100. 

The  influence  on  the  stress  peaks  due  to  change  in 
overlap  length,  2c,  is  similar  to  the  influence  of  E/Ea.  As 
c/tc  increases,  the  peak  stresses  and  stress  gradients 
decrease  in  magnitude. 

The  stress  singularity  parameters  together  with  the  order 
of  singularity  fully  describe  the  stress  fields  near  the  ends 
of  adhesive  overlap  and  at  the  bonded  corner  of  the  butted 
region  between  two  central  adherends: 

»ij/P«v  =  C„  (r/c)*-1 
where  c4J  are  the  stress  singularity  parameters  and  (p-i)  is 
the  singularity  order  (the  strength  of  the  singularity) . 
The    singularity   at    the   bonding    corners    of   two  different 
materals  depends  on  the  materal  constants  a  and  6  (Eguation 
5.16),   which  are  functions  of  Young's  moduli  and  Poisson's 
ratios  of  the  adhesive  and  adherends.      The  strength  of  the 
singularity   (   singularity  order)    a  =  p-i  decreases  as  the 
ratio,  E/E.,   increases.     The  singularity  order,  a,  increment 
is  very  small  over  the  range  E/Ea  >  50  and  is  noticeable  below 
the  range  E/Ea  <  30.     Efficient  double-doubler  joints  can  be 
designed  based  on  a  study  of  these  influences.  The  parameters 
of  the  joint  actually  designed  will  depend  on  the  particular 
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design  requirements  and  the  functions  the  joint  has  to 
accomplish. 

In  the  model  of  adhesive  bonded  joints  for  finite  element 
analysis,  it  is  important  that  the  critical  regions  of  the 
adhesive  be  modeled  with  very  small  elements.  It  has  been 
shown  that  the  stress  gradients  are  high  over  the  10%  of  the 
overlap  length  at  each  end  of  the  overlap.  Additionally,  5% 
of  the  overlap  length  at  each  end  of  the  overlap  and  15%  of 
the  adhesive  thickness  from  each  interface  must  be  modeled 
with  several  elements  in  order  to  get  more  accurate  peak 
stress  near  the  boundaries. 

The  specified  boundary  displacement  method  for  analyzing 
fine-mesh  submodels  in  the  critical  region  is  an  efficient 
method  for  solving  such  problem. 

6.3  Recommendations  for  Future  Work 

The   double-doubler   bonded   joint   with   linear  elastic 
adhesives  was  analyzed  in  this  study.    Future  study  should  be 
extended  to  joints  with  viscoelastic  adhesives  subjected  to 
constant  thermal  stresses,  transient  thermal  stresses  and  to 
dynamic   loads.       The   stress    singularity   parameters   may  be 
considered  as  a  means  of  failure  prediction  of  an  adhesive 
joint    by    comparing    the    calculated   parameters    with  known 
critical  values  of  these  parameters   for  the  material.  it 
should  be  pointed  out,  however,  that  this  kind  of  comparison 
can  be  made  only  for  machined  joints  where  contact  angles 
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between  the  adherends  and  the  adhesive  layer  do  not  vary 
(e.g., 90"  for  the  model  of  this  study). 

In  actual  manufacturing  processes,  however,  the  uniform 
contact  angles  do  not  occur  naturally.     Further  study  should 
be  done  to  evaluate  singularity  for  different  shapes  of  edge 
corner    with    identical    anisotropic    adherends.    Some  shape 
optimization     technigue     could     be     used     to     reduce  the 
singularity.     Furthermore,  the  stress  singularity  problem  is 
encountered  only  in  the  linear  elastic  stress  analysis.  If 
yielding  occurs  during  the  loading  process,  the  maximum  stress 
concentrations  will  be  significantly  reduced  due  to  plastic 
deformation.     This  should  be  further  studied. 
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Figure  A. 3 
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Figure  B.2  Case  1;  stress  ditributions  of  a  double-doubler 

joint  at  the  level  of  y=0.0675 
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Figure  B.4  Case  1.  stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.0925 
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Figure  C.l  Case  2.  Stress  distributions  of  double-doubler 
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Figure  c.2  Case  2.  Stress  distributions  of  double-doubler 

reinforcement:  of  a  continous  member  at  the  level 
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Figure  D.2  Case  2.  stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.0675 
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Figure  E.l  Case  3.  Stress  distributions  of  double-doubler 
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Figure  E.2  Case  3.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
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Figure  E.3  Case  3.  Stress  distributions  of  double-doubler 
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Figure  E.4  Case  3.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=-0 . 065 
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Figure  E.5  Case  3.  Stress  distributions  of  double-doubler 

of      OoS"1*  °f  3  COntinous  member  at  the  level 
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Figure  E.6  Case  3.  Stress  distributions  of  double-doubler 
reinforcement  of  a  continous  member  at  the  level 
of  y=0.07  J-eve± 
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Figure  E.7  Case  3.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 


I  T 

NORMAL  TzED  AX°IAL  COORDINATE''  X/C 


LU  D.03 

£ 
K 


|  ...»  . 

S 


LU   

< 

Q 

ft  *■ 


- 1  1  1  r 

0  *  0.  ft 


NORMALIZED  AXIAL  COORDINATE    x/c  ' 


1/3  o.i 

LU 
LT 


o.oe 


or 
< 

LU 

in 
a 

ft  - 


o 

2 


-r— t 
t.l 


i  l 
M 


151 


NORMALIZED  AXIAL  COORDINATE     X/C  ' 

Figure  E.8  Case  3.  Stress  distributions  of  double-doubler 
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of  y=0.095  61 
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Figure  F.l  Case  3;  Stress  distributions  of  a  double-doubler 
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Figure  F.2  Case  3.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=-0.08 
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Figure  F.3  Case  3.  Stress  distributes  of  a  double-doubler 

joint  at  the  level  of  y=-o.07 
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Figure  F.4  Case  3.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=-0.065 
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Figure  F.5  Case  3.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.065 
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Figure  F.6  Case  3.  Stress  distributions  of  a  double-doubl 

joint  at  the  level  of  y=0.07 
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Figure  F.7 
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Figure  F.8  Case  3.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.0925 
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Figure  G.l  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
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Figure  G.2  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=-0.08 
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Figure  G.3  Case  4.  stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=-0 . 07 
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Figure  G.4  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=-0. 065 
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Figure  G.5  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=0.065 
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Figure  G.6  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=0.07 


168 


D  CM 

co  o.o3 
id 

£ 


s 

s .... 

-o.o»  - 

HI 

I- 

Q 

111      -0  1 
IM 


NORMALIZED  AXIAL  COORDINATE  X/C 


I         Tj  1  1  1  1  ,  1  

NORMALIZED  AXIAL  COORD  I  NATE  '  'x/C 


T 

NORMALIZED  AXIAL  COORD  I  NATE "'x/C 


Figure  G.7  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=0.08 
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Figure  G.8  Case  4.  Stress  distributions  of  double-doubler 

reinforcement  of  a  continous  member  at  the  level 
of  y=0.095  1 
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Figure  H.l  Case  4.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=-o.095 
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Figure  H.2  Case  4;  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=-0.08 
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Figure  H.3  Case  4.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=-o.07 
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Figure  H.4  Case  4.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=-0.065 
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Figure  H.5  Case  4.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.065 
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Figure  H.6  Case  4.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.07 
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Figure  H.8  Case  4.  Stress  distributions  of  a  double-doubler 

joint  at  the  level  of  y=0.095 
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